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ABSTRACT 

Factorization properties of one-loop gauge theory amplitudes have been used as checks on explicitly 
computed amplitudes and in the construction of ansatze for higher-point ones. In massless theories, 
such as QCD at high energies, infrared divergences complicate factorization. Here we prove that 
factorization in such theories is described by a set of universal functions. In particular, a proof of 
the universality of one-loop splitting functions as the momenta of two particles become collinear 
is presented. Factorization in multi-particle channels is also given. The discontinuity functions 
that appear in the splitting functions may also be used to obtain infrared divergent box integral 
functions from finite ones. 
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1. Introduction 



In recent years a variety of one- loop gauge theory amplitudes with five or more external legs, 
including all virtual corrections to three-jet production at hadron colliders [1,2,3], have been com- 
puted. A number of formal developments in computational techniques, including spinor helicity 
methods [4], string-based methods [5,6], supersymmetry relations [7], recursive techniques [8,9], 
and unitarity methods [10,11] have allowed the computation of a large variety of new amplitudes. 

Factorization properties [12,1,13,10,3] provide strong consistency checks on newly computed 
amplitudes. Furthermore, they provide a means for obtaining higher-point amplitudes from lower- 
point ones. To do this one constructs a function that has the correct factorization properties in all 
channels, which is then an ansatz for the amplitude. An explicit construction of n-gluon amplitudes 
with all identical helicities derived from this technique has been given in refs. [13,14] and verified 
via recursive techniques [9]. 

The limiting form of one-loop massless gauge theory amplitudes as the momenta of two legs 
become collinear has been given, although not proven, in terms of splitting functions in previous 
papers [13,15,10,3]. In particular, a tabulation of the splitting functions in QCD has been given 
in refs. [10,3] from the collinear limits of calculated five-point amplitudes. Ample evidence of the 
universality of the splitting functions is provided by the large number of explicitly known one-loop 
helicity amplitudes in massless gauge theory [5,16,1,2,14,13,9,10,11]. For multi-particle factorization 
less information is available; except for the recently computed set of one-loop six-point amplitudes 
in N = 4 supersymmetric gauge theory [11], no examples have been constructed with non-trivial 
factorization properties. (Mahlon [9] has also constructed a set of helicity amplitudes containing 
multi-particle poles, but with simple factorization properties due to a lack of infrared singularities.) 
Other properties at loop level which are known are for the unpolarized differential cross-sections, and 
give the well known Altarelli-Parisi splitting functions [17]. The Altarelli-Parisi splitting functions 
are for cross-sections which include both real emission and virtual diagrams; the splitting functions 
that we discuss in this paper are for the virtual (i.e. loop) diagrams alone. 

In this paper we provide a general proof of the universal behavior of massless amplitudes, in- 
dependent of the number of legs, when factorized on a collinear or multi-particle pole. Amplitudes 
which lack infrared divergences factorize straightforwardly as one would expect. For amplitudes 
with infrared singularities the situation is more subtle: massless gauge theory amplitudes generally 
do not factorize in any simple sense since there are contributions which cannot be interpreted di- 
rectly in terms of lower-point amplitudes. Nevertheless, we will prove here that the non-factorization 
may be described in terms of a small set of factorization functions, given in this paper, whose coeffi- 
cients are fixed by the known infrared divergences. We also give a practical procedure for calculating 
the splitting functions of any one-loop amplitude from two- and three-point diagrams. 

In previous papers [13,15] proofs of the universality of the (two-particle) one-loop splitting 
functions for the case of n external gluons with a virtual scalar loop were outlined. (This was suf- 
ficient for the particular helicity amplitudes considered in these papers, because of supersymmetry 
identities [7] which relate the contributions of virtual gluons and fermions to those of scalars.) The 
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proof was based on an analysis of Feynman diagrams in the limit that the momenta of two external 
legs become collinear. This type of proof also works for a fermion loop [18], but a generalization 
to the case of virtual gluons and multi-particle factorization is non-trivial due to complications in 
quantifying non- factorizing contributions. 

One expects, however, that the universality of the splitting functions should follow from general 
field theory considerations and not from the details of particular diagrams. Our proof will make this 
explicit and is valid for any particle content and external states. The key behind the proof is that we 
can link all non-factorizing contributions to infrared divergences present in the one-loop amplitudes. 
We will show that only a small number of functions may enter and that their coefficients are fixed 
by their infrared divergences. Since the divergences have a known form [19,20,21], all functions 
entering into factorization are fixed and universal; this is the main result of this paper. 

A spin-off from our analysis is that the same functions that contribute to the loop splitting 
functions may also be used to efficiently obtain all infrared divergent box integrals by taking limits 
of known infrared finite massive box integrals [22]. To illustrate this, we reproduce the box integrals 
necessary for calculating amplitudes in massless gauge theory given in ref. [23]. 

First, in Section 2 we review some properties of gauge theory amplitudes that we will find 
useful in our subsequent discussion. In Section 3 we present our results for the one-loop splitting 
functions; an explicit sample calculation is presented for obtaining the splitting functions from 
three-point diagrams. In section 4 we present analogous results for multi-particle factorization 
and a non-trivial example. The actual proof of these results are presented in Section 5 and the 
appendices. In section 6 we show that the same functions which appear in the non-factorizing 
contributions are useful for obtaining infrared divergent box integral functions from infrared finite 
ones. We also provide a variety of appendices on integrals and their properties which are useful in 
the proof. 

2. Review of Previous Results 

We now briefly review some previous results and conventions that we will use in this paper: 
the color ordering of non-abelian scattering amplitudes, the spinor helicity method, tree-level two- 
particle collinear limits and the structure of infrared and ultraviolet singularities. We review (and 
slightly modify) the integral reduction method for calculating any one-loop integral in terms of box, 
triangle and bubble functions in appendix I. 

2.1 Color-Ordered Amplitudes 

Tree-level SU(N C ) gauge theory amplitudes can be written in terms of independent color- 
ordered partial amplitudes multiplied by an associated color trace [24]. This is extensively discussed 
in the review article of Mangano and Parke [25], whose normalizations and conventions we follow. 
(In particular, we normalize fundamental representation color matrices as Tr[T a T b ] = 5 ab .) One of 
the key features of the partial amplitudes is that the external legs have a fixed ordering. 

At one- loop, although the analogous decomposition [26] is a bit more complicated, gauge 
theory amplitudes may be conveniently written in terms of gauge-invariant 'primitive' amplitudes 
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[3] which also have fixed ordering of external legs. The n-gluon primitive amplitudes correspond to 
the leading-color partial amplitudes, ^4^i P . Furthermore, the complete amplitude may be expressed 
in terms of appropriate permutation sums over primitive amplitudes multiplied by ordered color 
factors. Therefore, all results obtained in this paper for primitive amplitudes may be converted to 
results for full amplitudes. For amplitudes with only external gluons, or two quarks and the rest 
gluons, explicit expressions giving the full amplitudes in terms of these primitive amplitudes may be 
found in refs. [10,3]. For this paper, the main property needed is that a decomposition of one-loop 
SU(N C ) gauge theory amplitudes exists in terms of a set of color decomposed amplitudes where 
the ordering of external legs is fixed; this is convenient, although not necessary for our discussion. 

2.2 Spinor Helicity 

In explicit calculations it is usually convenient to use a spinor helicity basis [4], where all 
quantities are rewritten in terms of Weyl spinors |/c ± ). In the formulation of Xu, Zhang and Chang 
the polarization vectors are expressed as 

V ' V2(qk) V > y/2[kq] V ' 

where q is an arbitrary null 'reference momentum' which drops out of the final gauge-invariant 
amplitudes. The reader is referred to the article of Xu, Zhang and Chang for further details. For 
the purposes of this presentation we note that 

(K\k+) = (ij) = y/2k t • k, exp(^), {kf\kj) = [ij] = -v^fci-fyexpH^) , (2.2) 

where is a phase. These spinor products vanish in the limit fc, • kj — > 0, are anti-symmetric and 
satisfy 

(ij)\ji]=2k i -k j . (2.3) 

2.3 Tree-Level Two-particle Collinear Factorization 

Consider first an n-point tree-level partial amplitude A^ ee (l, 2, . . . , n) with a fixed ordering of 
external legs and an arbitrary helicity configuration (Ai, A2, • • • , A n ), where each Aj = ±. As the 
momenta of two neighboring legs a and b become collinear the leading behavior of the amplitudes 
is given by [27,12] 

^tree ^ ^ Split^ (fl A » , b Xb ) A^\{. . . K X . . .) , (2.4) 

A=± 

where the non- vanishing splitting functions (or amplitudes) diverge as 1/ ^/s^ in the collinear limit 
s a b = (k a + kb) 2 — ► 0. In this equation and all subsequent ones we extract all coupling constants 
from the color-ordered amplitudes. The collinear limit is defined by k a = zK and k^ = (1 — z) K, 
where the null vector K is the sum of the collinear momenta; A is the helicity of the intermediate 
state with momentum K. The tree splitting functions Split !!^ e (a Aa , b Xb ) may be found in ref. [25]. 
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The g — ► gg splitting functions are 



Split!! ee (a" , b~ ) = , Split!! 66 (a+ , 6+) 



Split!! 66 (a+ , 6" ) = ^ , Split!! 66 (a" , b + ) - ( J " : ' " 



(2.5) 



V^(l-z)[a6] ' : ' " ' Vz(l-z)[a6] ' 



which explicitly exhibit the \/y/s a b pole. All remaining g ^ gg splitting functions may be obtained 
by parity. 

2.4 Infrared and Ultraviolet Singularities in One-Loop Amplitudes 

In massless gauge theory one encounters three types of singularities when evaluating virtual 
corrections with fixed particle number. Besides the usual ultraviolet singularities one encounters 
both soft and collinear infrared ones. (By 'soft' we refer to the e -2 terms in dimensionally regulated 
amplitudes.) These types of singularities must cancel in final physical cross-sections containing both 
real and virtual corrections, but are present in the individual parts. The singularities in an n-point 
massless QCD amplitude with a fixed ordering of legs are of the form [19,20,21] 



^loop 



singular 



a tree 



i n / 2 \ e 1 



(2.6) 



where 



(4vr) e r(l + 6)r 2 (l-e) 
Cr " 16vr2 r(l - 2e) ' (2 - 7) 



with e = (4 — D)/2 the dimensional regularization parameter and fi the renormalization scale. The 
parameters tSj™' are the coefficients of the soft singularities and is the sum of the coefficients 
of collinear and ultraviolet divergences, which depend on the particle content of the amplitude. 
Although we will write the soft singularities in the form appearing in eq. (2.6), the functional form 
of the amplitudes are taken, as usual, to be expanded in e. This expansion is implied in eq. (2.6) 
and in all subsequent expressions. 

The main property that we will make use of is that in massless QCD the coefficients <sj n ' and 
C' n ' are known. For primitive amplitudes which have been stripped of all color factors and have a 

[nl 

fixed ordering of legs, these coefficients are particularly simple. For such amplitudes, Sj is or 1 
depending on the particle type of the internal loop line connecting legs j and j + 1 in the 'parent 
diagram', which is the one where all external legs are directly attached to the loop by three- vertices. 
If the propagator between legs j and j + 1 is a gluon then = 1, while if it is a fermion or scalar 
then tSj"^ = 0. (Explicit examples are found in ref. [3].) 

The collinear infrared singularities for n-point amplitudes are 

n 
a=l 
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where the sum is over all n legs and depends on whether leg a is a gluon g or a fermion q, 

11 1 rif 1 n s 



6iV c ' 



(2.9) 



7(f) 

^)-!('-^)- 

where rif is the number of fermions, n s the number of scalars (which is zero in QCD) and N c the 
number of colors. Our conventions follow the ones in refs. [1,10,3]; in particular we have extracted 
an overall factor of iV c from leading color amplitudes so that contributions from fundamental 
representation loops carry a factor of 1/N C . The normalization for each complex scalar is non- 
standard and represents a total of four states in the (N c + iV c ) representation (instead of the usual 
two; the rationale for this choice is to maintain supersymmetry identities with the four states of 
fundamental representation Dirac fermions). There are also ultraviolet singularities 

11 I rif ln s 



C H -(n-2)(---^--^) 
The coefficient appearing in eq. (2.6) is the sum 

An] _ An] , An] 
u ~~ U UV U IR ) 

of the infrared and ultraviolet contributions in eqs. (2.8) and (2.10). 



(2.10) 



(2.11) 



3. Collinear Limits 

3.1 General Results 

In this section we present our results for the factorization of one-loop amplitudes as the mo- 
menta of two color-adjacent external legs become collinear; in subsequent sections we present the 
proof. The behavior of one-loop amplitudes in the collinear limit found from explicit calculations 
of amplitudes with five or more legs, is [1,10,3] 

A loo P ^ y- I g plit tree (fl A ^ ft A b ) ^loop ^ _ _ R X _ _ ^ + g^op, . A 



y: {spht* 

A=± ^ 



(a 



(3.1) 



where the A l ° op and A^ ee are color-decomposed one-loop and tree amplitudes with a fixed ordering 
of legs and a and b are consecutive in the ordering. This is a natural generalization of the tree-level 
factorization and is exactly what one would expect in massive theories with no infrared divergences, 
since in this case individual diagrams respect the factorization depicted in fig. 1, where legs i and 
i + 1 correspond to a and b. 





+ 




Figure 1. The 'naive' factorization of one- loop amplitudes in the limit as two external 
momenta become collinear, where K = ki + fcj+i. The shaded disc represents the sum over tree 
diagrams and the annulus the sum over one-loop diagrams. 
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For infrared divergent theories the situation is much more subtle since individual diagrams may 
not have a smooth behavior as the intermediate momentum becomes massless. We shall, however, 
prove that the splitting functions Split are independent of the number of external legs n, and 
that they may be obtained by calculating three-point diagrams followed by adding in a particular 
set of functions fixed by the singularities in e through eq. (2.6). 

The following analysis is performed on bare amplitudes (before renormalization) . For massless 
amplitudes it is slightly more convenient to leave all quantities as bare, but at any point in the 
discussion one may convert to renormalized quantities by performing the appropriate ultraviolet 
subtraction in all quantities; this will be done for the final splitting functions. We will also use 
the dimensional reduction [28] or equivalently the four-dimensional helicity (FDH) regularization 
schemes [5]. These schemes are convenient since they do not introduce any e-dependence coming 
from the contraction of tensors which count the number of states circulating in the loop. It is 
straightforward to convert the splitting functions to conventional [29] or 't Hooft-Veltman schemes 
[30]; through 0{e ) the splitting functions shift by constants multiplied by tree splitting functions, 
as given in refs. [10,3]. A discussion of scheme conversions has been given in ref. [16]. 

In general, infrared divergent amplitudes will not factorize naively. As a simple example of 
non-factorization, consider an infrared singularity of the form 



1 



€ 



1 1 



c r ^(-s 23 )^AT e = ~c r - j - - ln(- S23 ) + • • • < ree , (3.2) 



e 2 e 



in the collinear limit k\ = zK and k 2 = (1 — z)K. This singularity contains the logarithmic term 

l\n(-s 23 )^^\n(-(l-z)s K3 ), (3.3) 

with kf = 0, which introduces a ln(l — z)/e not belonging with either the 'naively factorized' 
diagrams on the left- or right-hand-side of the tree-pole in fig. 1 (for i = 1). It therefore cannot be 
interpreted as a factorizing contribution. 

Our results for the limits of amplitudes as legs i and i + 1 become collinear (ki — > zK and 
k i+ i — > (1 — z)K) are as follows. As the two legs become collinear, loop integrals may not have 
smooth limits and develop infrared divergences. The possible discontinuities, which describe the 
off-shell to on-shell transition, are described by a set of universal 'discontinuity functions'. Two 
discontinuity functions which play an explicit role in computing the loop splitting functions are 

u \ 1 ( V 2 



f7^y 1 ' (3 ' 4> 

(There are additional discontinuity functions to be discussed in section 5, but they are not needed 
for our discussion of the factorizing diagrams.) 

Now consider the computation of loop splitting functions which are composed of factorizing 
and non-f actor izing pieces, 

Split loop = Split fact + Split non - fact . (3.5) 
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The first step in obtaining the factorizing contributions is to compute the diagrams depicted 
in fig. 2. Observe that diagrams with bubbles on external lines are not included. This is due to the 
dimensional regularization prescription that massless on-shell bubble diagrams vanish [31], which 
is interpreted as a complete cancellation of infrared and ultraviolet divergences. (Below we also 
describe the diagrams necessary for calculating the splitting functions in a massive theory.) The 
diagrams in fig. 2 yield the general form 



V = Bo 



1 



Si,i+1 



+ 5!- + So, 

e 



(3.6) 



where Bi and B 2 are rational functions depending on the particle content and on the type of 
external legs. Depending on the type of off-shell leg, represented by the dotted line in fig. 2, V may 
have uncontracted spinor or vector indices. 

i + 1 



+ 




>o~ + >o~ 



Figure 2. The diagrams in a massless theory (ignoring tadpoles) that need to be calculated 
to obtain the factorizing contribution to the loop splitting function. The dotted line represents 
the off-shell leg on which the collinear factorization is performed. 

It is convenient to collect all singular terms, including those in eq. (3.6), into the non- factorizing 
category. We achieve this by subtracting all singularities in e from eq. (3.6) using the discontinuity 
functions (3.4), and then adding them back into the non- factorizing set discussed below. This 
collects all non-smooth behavior into the non-factorizing contributions. (In general, the coefficients 
of the singularities, B\ and B 2 , depend on gauge choices and the processes under consideration.) 
Subtracting out the divergences from the factorizing diagrams yields 



V' = Bo 



A* 



i+i, 



+ B X 



- - b(si, i+ i) 



+ B 



(3.7) 



which we label as the 'factorizing' contribution to the splitting functions. By construction this is 
completely free of singularities in e. The result V' may then be contracted against the (n — 1)- 
point tree diagrams, as in the last diagram in fig. 1, and converted to an expression in terms of 
spinor helicity. The conversion is performed by inserting a complete set of helicity states on the 
intermediate factorized leg and then taking the collinear limit; this yields 



^Sph^M + iK^U--^---) 



(3.8) 



In this way we obtain the factorizing contribution to the loop splitting functions. In the next 
sub-section we present an explicit example. 
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The second type of contribution is the 'non-factorizing' one, which comes from non-smooth 
behavior in any of the diagrams. This includes the non-smooth contributions which were subtracted 
from the diagrams in fig. 2 and any other diagrams which generate kinematic poles from the loop 
integrals. These contributions are proportional to discontinuity functions (3.4) and a limited set of 
integral functions containing poles in s^j+i to be discussed in section 5. 

We find that the non-factorizing contributions to the loop splitting functions are proportional 
to the tree- level splitting functions, so we have 



Split™T fact (a\k Ab ) 



c r x Splits (a x «,b Xb ) x r s (a,b) , 



(3.9) 



where r s contains no spinor products or helicity dependence, and a and b are legs i and i + 1. 
The contributions to r s are given in Table 1 in terms of the singularities of the amplitudes (2.6). 
The proportionality of Split non_fact to Split tree follows from the appearance of singularities in all 
non-factorizing contributions and from the proportionality of the singularities in loop amplitudes 
to tree amplitudes. 

In Table 1 the coefficients for the singularities <Sj"i, S\ n \ and are for the n-point 

amplitude, and the coefficient Cl™ -1 ! is for the (n — l)-point loop amplitude described by the second 
set of diagrams on the right-hand-side of fig. 1. The additive contribution to r s is given by the 
coefficient in the first column multiplied by the corresponding terms in the third column. Note that 
r s depends only on the particle types of legs i, i + 1 (and the fused leg) through the dependence 
onS\ n] andCM - & n ~ l \ 



Coefficient 



Singularity 



Non-Factorizing Contribution to r s 



-Li 2 (l-z) 



S. 



i+l 



H+1A+2 



-(i-«K< 



+i 



Li 2 (z) 



-4(^) 



-^(— ) 



1 



e(l-2e) V-s 



2,2 + 1 



Table 1: The 'non-factorizing' contributions to the loop factorization functions in the S^j+i 

channel. The coefficients in the first column are the coefficients of the contributions to r s given 
in the third column. 

The first two entries appearing in the third column are the collinear limits of discontinuity 
functions summed with specific box functions 
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fcj||fci+i 

ki\\k i+1 



(3.10) 



where F^ +2 and F*™ +3 are the single-external-mass (or off-shell) box functions given through 
O(e ) in eqs. (IV.14e) and (IV.lOe). (The form in eq. (3.10) is valid to any order in e.) 

The terms of higher order in e can become important when performing phase-space integrals 
for n + 1 parton contributions to n-jet final states at next-to-next-to-leading-order (NNLO) when 
using, for example, the formalism of ref. [20]. The splitting functions may be used to compute 
analytically the phase-space integrals in the collinear regions; in these regions the phase space 
integrals contain powers of e _1 which can cancel against higher order terms in the splitting functions 
to leave finite results. The higher-order in e terms in the splitting functions may be obtained (within 
the dimensional reduction [28] or FDH [5] schemes) by simply keeping as many orders as desired 
from both the diagrams in eq. (3.7), and from Table 1, using the form in eq. (3.10) for the first 
two entries. The other functions appearing in the table are given in a form valid to all orders in 
e. In performing the phase space integral it is convenient not to convert to a helicity basis, but 
to leave all expressions in terms of the formal polarization vectors. (If one were to convert to a 
helicity form one would have to account for the fact that the momentum over which one performs 
the phase space integral is actually in 4 — 2e dimensions and not in four dimensions; this type of 
subtlety has, however, already been addressed by Mahlon [9] in the context of recursion relations 
for amplitudes.) 

In summary, the total contribution to the loop splitting function (3.5) is given by the sum of the 
factorizing and non-factorizing contributions. The factorizing contributions (3.8) are independent 
of the number of external legs since they always come from the same three-point diagrams depicted 
in fig. 2. The non-factorizing contributions (3.9) are also universal because their coefficients are 
fixed by the singular terms (2.6). Thus, for any number of external legs the loop splitting functions 
appearing in eq. (3.1) depend only on the 'local' properties within the diagrams: the helicity and 
particle type of the two collinear legs and the virtual matter content. For the color-ordered primitive 
amplitudes of ref. [3] the splitting functions also depend on the routing of the fermion through the 
diagrams. 

3.2 Collinear Limit Example 

As an example, consider a one-loop n-gluon partial amplitude in a theory with rif massless 
fermions, n s massless complex scalars, and iV c colors in the limit that legs 1 and 2 become collinear: 
k\ = zK and k 2 = (1 — z)K. First we compute the factorizing contributions by evaluating the 
two- and three-point loop diagrams of the type in fig. 2 (for i = 1); the ones with fermion loops are 
depicted in fig. 3. The gluon and scalar loop diagrams are similar except that one must include 
diagrams with four-point contact vertices. In performing the calculation, the intermediate leg (the 
dashed line) with momentum (— k\ — k 2 ) should be left off-shell in the initial part of the calculation. 
These diagrams are conveniently evaluated [32] using color-ordered Feynman background field gauge 
[33] for the loop and Gervais-Neveu gauge [34] for the tree parts of the diagrams. Through all orders 
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in e we have 



Background 



i T r / n s n f \- 



ei-A 



2 t-2 



fei • ko 



where 



r r = 



6 r(l + e)r 2 (l-e) 



16vr 2 



+ 0(e) 



(3.11) 



(3.12) 



(4vr) 2 - e T(4-2e) \s 12/ 

In the computation we have used the dimensional reduction scheme, while in the conventional [29] 
or 't Hooft-Veltman [30] schemes there would be an extra overall factor of (1 — e) in the gluon loop 
contributions, which has no effect through O(e ). The fermions and scalars are taken to be in the 
fundamental representation. (For the adjoint representation the n s and n/ terms would not have 
a factor of 1/N C .) A background field Ward identity between the two and three-point function 
cancels the divergences between the separate diagrams leaving the finite result in eq. (3.11). 

One may of course use other gauges, although in general these yield more complicated expres- 
sions. For example, with color-ordered Feynman gauge the gluon loop diagrams are 



Feynman yf^, I 1 [ /J. 



+ i c r b(s 12 ) -ei ■ e 2 (k 1 - k 2 y + ■ e 2 (h + fc 2 ) M + e 2 ■ k^ - e 1 ■ k 2 e% 



+ 3 L 



£1 • £2 



£l • h 



h ■ k 2 



i-k 2 y + -(3-2e)(k 1+ k 2 y 



(3.13) 

which has been written in a form to expose the discontinuity functions (3.4) contained within 
the diagrams. The result for this gauge (3.13) looks rather different than the one for background 
field gauge (3.11) and is a reflection of the gauge dependence of the diagrams due to the off-shell 
intermediate leg. 





+ 
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Figure 3. The diagrams for the factorizing massless fermion loop contributions to the 
gg loop splitting functions. 



As discussed above, the first step in obtaining the contributions to the loop splitting function 
after having calculated the factorizing diagrams in fig. 2 (or fig. 3 for fermions) is to push all singu- 
larities in e into the non-factorizing category using the discontinuity functions (3.4). After moving 
the singularities from the Feynman gauge result (3.13) into the non-factorizing contributions, it 
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agrees with the background field gauge result (3.11), up to terms proportional to (ki + k 2 ) fl . Since 
the free index fi contracts against a tree- level conserved current, these terms are irrelevant. Thus, 
in either gauge, we obtain the 'factorizing' contribution as 



.fact W (, , n s n f \ t t 
n \ + W c W c )~6 

where 



£1 • k 2 e 2 ■ ki 

£l ■ £ 2 - 



{k!-k 2 r_ 4tree 



r lliV A^ 1 {K,S,...,n) v , (3.14) 



2 



^n-i(^ 3, . . . , nf = -^L^A* e _\(K, 3, . . . , n) . (3.15) 

and i]^ is the Minkowski metric. Although we obtain the same result in either gauge, the natural 
choice for the calculation is background field gauge, since one directly obtains the physical part from 
eq. (3.11). (The form in eq. (3.14) is valid to higher order in e and thus would be the appropriate 
expression to use for NNLO phase-space integrals.) 

The expression (3.14) may be converted to spinor helicity notation by inserting a complete 
set of helicity states on the intermediate leg [8] using an off-shell generalization of the helicity 
formalism. For on-shell momentum K from eq. (2.1) we have 

(3 . 16 ) 

2(<? \fC\q ) 

The right-hand-side of this equation is well defined for off-shell K, so we may take it as the definition 
for a product of off-shell polarization vectors; we need only define the product of polarizations since 
this is the only combination that appears when inserting a complete set of helicity states. It is not 
difficult to verify with eq. (3.16) the identity 

= -e+(-K, q) e- (K, q) - e~ (-K, q) £+ (K, q) + K ^ + ^ K - , (3 .17) 

A • q 

which then allows us to replace the r/^ with off-shell polarization vectors plus a piece which will 
vanish by current conservation. For K = k\ + k 2 we have 



(qk 2 ) [k 2 h] (ki 






- (qh) [ki k 2 ] (k 2 






2(q- 







(k 2 - h) ■ e + (-K, q) e~ (K, q) = " L * " w 1 L ^ 1 1 . (3.18) 

We parametrize the collinear limit by k\ — ► zK and k 2 — ► (1 — z)K with K 2 — ► 0, to obtain 

(fc - fci) • e + (-K, q) e- (K, q) - -^(l-z) [A* fc 2 ] e; (Jf, ?) , (3.19) 

where the polarization vector on the right-hand side is on-shell. Thus after inserting a complete 
set of helicity states using eq. (3.17) (through C(e )) we obtain from eq. (3.14) 

A f : ct ^ { 1 + jr~¥) E^^ 2 ^™^ 3 n )> ( 3 - 2 °) 

c c X=± 
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where 



Split f ^ ct (l+,2 



Split f ^ ct (l",2+; 



0, 



S P lit f r(l + ,2+) = -^y/^z)^ , 



(3.21) 



Split f _ act (l + ,2 + ^ 



1 



(12)' 



are the factorizing contributions to the loop splitting functions. 

To obtain the non-factorizing contribution to the splitting function we note that for n external 
gluons (S] n ' = 1 for all i and = C^™ -1 ], so that from eq. (2.6) the singular terms in e are 



A loop 



singular 



a tree 



1 ^ / ii 2 V 2/11 In/ 1 n s \ 



(3.22) 



From Table 1, the non-factorizing contribution to r s for this example is the sum of the first three 
entries in the last column, given by 



r s (l 



1 



-ZSV2 



1 ( ( /\ Y - Li 2 (l -z)- Li 2 (z) + 0( C ) . (3.23) 
e 2 V-(l - z)s 12 / 



After simplifying r s using the dilogarithm identity 



7T 



Li 2 (l - z) + Li 2 («) = - ln(z) ln(l - *<) + — , 
we obtain from eq. (3.9) the non-factorizing contributions to the g ^ gg splitting functions, 



(3.24) 



Split 



non— fact 



(l,2)=c r Spht t _T(l,2) 



1 



7T 



+ 21n(z)ln(l -z)- — 



e 2 y z (l - z)(-s 12 ) 
Combining the results (3.20) and (3.25) yields the total contribution 

Split loop (l, 2) = (l + ^ - ^t)s P lit fact (l, 2) + Split non - fact (l, 2) , 



(3.25) 



(3.26) 



in agreement with the results in ref. [10] obtained from taking the explicit collinear limits of five- 
gluon amplitudes [1]. 

Finally by performing the ultra-violet subtraction using (2.10) we obtain the corresponding 
splitting function for renormalized amplitudes as 



Splitj£ p (l,2) = (l + ^ - ^)Spht fact (l,2) + Spht non - fact (l,2) 

-c r I(H_IZV _IZM S plit tree (l,2), 
r e V 6 3iV c 6N C J 1 y ' ; ' 



(3.27) 



where all remaining powers of e _1 are infrared singularities. 

The loop splitting functions for amplitudes with external fermions and gluons may similarly be 
obtained; these have already been extracted from four- and five-parton calculations and tabulated 
in refs. [10,3]. 
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3.3 Collinear Factorization with Massive Fermions and Scalars 

For amplitudes with no infrared singularities the splitting functions are entirely determined by 
the naively factorizing contributions. Consider, for example, the collinear limit of gluon amplitudes 
in a theory with massive fermion or scalar loops. Since the mass cuts off any singular behavior 
as legs go on-shell, all loop integrals are smooth in the collinear limit and the amplitude naively 
factorizes; the collinear poles only come from external gluon propagators. Thus massive fermion 
contributions to the splitting functions may be directly determined by calculating only the diagrams 
in fig. 4. (For a scalar loop there are additional diagrams containing the four-point interaction.) 

Following the usual renormalization procedure, half of each bubble on the external line is asso- 
ciated with external wavefunction renormalization. Note that the set of diagrams that contribute 
to the massive loops is different from that of massless loops (given in fig. 3). The difference between 
the two sets is due to the dimensional regularization prescription that massless bubbles on external 
legs vanish; in fig. 3 the two diagrams with bubbles on external lines vanish and are therefore 
not included. For the massive case in fig. 4 half of the bubble on the 'internal' line on which the 
factorization is being performed belongs with the (n — l)-point amplitude and not with the loop 
splitting function. 



« + l 




Figure 4. The diagrams representing contributions of massive fermions to the g — > gg loop 
splitting function; the wavy lines represent gluons, and the dashed lines represent the fact that 
only half of each bubble is included. 

The calculation of the diagrams in fig. 4 is straightforward and gives the entire contribution 
to the loop splitting functions. This gives the contribution of massive fermion and scalar loops to 
the g — > gg splitting functions as* 



Spht!_°°>,&) = - 




where n/ and n s are the number of fermions and scalars and m/ and m s the corresponding masses. 

The proportionality of these contributions to the tree splitting functions may be understood 
from the structure of the effective action (for k\ ■ k 2 <C m 2 ). The contribution to the effective 
Lagrangian from a massive fermion or scalar is 

C eS = -^F^F^ + 0(-^), (3.29) 

* As in refs. [1,3], each scalar here contains four states (to match the four states of Dirac fermions) so that n s 
must be divided by two for comparisons to conventional normalizations of scalars. 
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where is the wavefunction renormalization. The higher order terms are suppressed in the 
collinear limit because they contain l/m 2 instead of 1/fci • kz, leaving only a renormalization of the 
tree splitting functions. In the massless case this argument breaks down and there is no reason to 
expect the loop splitting functions to be proportional to the tree ones. 

Subtracting off ultra-violet singularities (2.10) from eq. (3.28) gives the splitting functions from 
massive fermions or scalars for renormalized gluon amplitudes as 



Spht_ A ( a ,b) = -- 



n f 



In 



3 (4tt) 2 [N c \m 2 J 2N ( 



Split t _ r r(a,6) + C(e) 



(3.30) 



The gluon loop contributions in this theory may be obtained from eq. (3.27) by taking the number 
of massless scalar and fermions to vanish, nj = n s = 0. 



4. Multi-particle Factorization 

Consider now multi-particle factorization which we will show has analogous behavior to the 
two-particle collinear limits. As in the latter case, this non-smooth behavior may be linked to the 
infrared singularities appearing in massless amplitudes. 



4.1 General Considerations 

We will prove that the factorization properties for (A^ + k i+ i + 
(with r > 2) are described by the universal formula, 



+ ^i+r-l) 2 



K 2 



j ki 



+^r+l(^«5 ' ' ' ki+i — 1, K ) 2 A n _^, +1 (K , /Cj_|_ r , ■ ■ ■ , ki-i) 

-\-A r _ Sr i{ki 1 ■ ■ ■ ki+ r —i, K ) A n _ rJrl (K ,ki+ r ,---,ki-i)c-pJ- n (K ;ki,...,k n ) 

(4.1) 

where the one- loop factorization function T n is independent of helicities. This formula is similar to 
the one for an amplitude which factorizes naively, as depicted in fig. 5, except that T n may contain 
kinematic invariants with momenta from both sides of the pole in K 2 ; for example ln(— t^j) = 
ln(— Si-i,i) is one such logarithm. (In fig. 5 we have made the bubble on the intermediate leg explicit; 
in fig. 1, for collinear limits, the bubble diagram was implicitly included in the last diagram on the 
right-hand-side, as depicted in fig. 2.) As for the splitting functions the factorization function is 
composed of factorizing and non-factorizing components, 



.r-fact _|_ r-non- fact 



(4.2) 



For convenience we have extracted an overall factor of c r from the factorization function. 
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i + r - 1 




Figure 5. The diagrams with an explicit pole in K 2 = tf^ coming from a tree propagator. 

One of the results that we shall prove is that only a few integral functions may contribute to a 
pole in any multi-particle channel. Furthermore, the singular terms tightly constrain the coefficients 
of the allowed integral functions. The factorizing contributions in a multi-particle channel are 
composed only of integral functions depicted in fig. 6; all external momenta from one of the two 

\r] 

sides of the multi-particle pole are part of one leg of these integral functions. If t\ constitutes a 
kinematic invariant of an entire leg of an integral function there will be non-factorization from the 
discontinuous limit. There are also two integral functions which may enter into non-factorization: 
the box functions i4 :7 — i;i+i and Lt :ra _ r _i ; j_|_ r _|_i in fig. 7. (The values of these box functions are 
given in appendix IV.) These two box functions have the property that they contain explicit poles 
in K 2 yet they would not appear in any of the naively factorized diagrams in fig. 5. 




: m 




Figure 6. The integral functions which contribute to diagrams with a tree pole in the 
K 2 = tP channel. 



■; + t - 1 



: + r-l 



=7 ' V 



(a) (b) 

Figure 7. The two boxes functions which contribute to non-factorization in the channel t\ , 
indicated by the dashed line. 

The procedure for obtaining the factorization function is similar to the one for obtaining the 
loop splitting functions. First compute the factorizing contribution given by the bubble loop in the 
third class of diagrams on the right-hand-side of figure 5. This contribution is of the form 



Ci-+C . 

e 



(4.3) 



16 



In this case there is no e~ 2 singularity since bubble diagrams have at worst e _1 ultraviolet singu- 
larities (in a Feynman-like gauge). Note that the bubble loop is gauge dependent. As was the case 
for two-particle factorization we must subtract an appropriate discontinuity function to obtain 

+ C , (4.4) 

which is finite. After inserting a complete set of helicity states we obtain the factorizing contribution 
jjrfact f romeq . (4.4). (i n 

some cases, such as factorizing on a fermion line, this contribution vanishes.) 
The second type of contribution to the factorization function, ^"non-fact ^ are ^ e non .f ac torizing 
ones given in Table 2. The coefficients of the contributing terms in the third column are given in 
the first column. As for two-particle collinear limits the non-factorizing terms are fixed by the 
singular terms (2.6). The box functions F 2mh appearing in the third column of Table 2 are defined 
in eqs. (IV. 10) and (IV. 14) and correspond to the kinematic configurations depicted in fig. 7. 



Coefficient 


Singularity 


Factorization Function Contribution 


cH 


i/ ^ y 


2(M 2 ) £ (« 1;i+1 + ^(-4 r] )- £ ) 


cH 

°i+r-l 


- U " 2 )' 

e Si-\-i — l,i+r ' 


2(/i 2 ) e ^F 2 ™^ r _ 1;i+r _ | _ 1 + — (— tj ') 




1 

e 


^ > e (i _ 2e) [ h > 



Table 2: The 'non-factorizing' contributions to the factorization functions in the channel t\ — > 
0. The three coefficients in the first column are the coefficients of the contributions to the 
factorization function given in the third column. 

The total contribution to the factorization function T n is the sum of terms in Table 3 along 
with ones in eq. (4.4). 

4.2 Multi-Particle Factorization Example 

We now present an example to illustrate the use of multi-particle factorization to provide 
checks on explicitly calculated results. This is analogous to the types of checks that have been 
performed using collinear limits both for tree [27, 12] and loop [13,10,3] amplitudes. Multi-particle 
factorization can also be helpful for constructing ansatze for higher-point amplitudes. 

Consider the one-loop N = 4 supersymmetric six-gluon amplitude A^ 1 =4 (l + , 2 + , 3 + , 4~ , 5~, 6~), 
which has been computed in ref. [11]. In the = (ki + k% + k^) 2 channel, the right-hand-side 
of eq. (4.1) vanishes since A tiee (1 + , 2 + , 3 + , K =t ) = 0. Thus the amplitude does not contain a pole 

f3l f3l 

in t\ . The t 2 channel, however, will contain a pole, since tree and loop amplitudes appearing 
on the right-hand-side of eq. (4.1) do not vanish. Since we are dealing with = 4 super- Yang- 
Mills, only box functions may enter [10] in the integral reduction. Therefore, the only integral 

[31 

functions which may contribute to the coefficient of the t 2 pole are given by the four box functions 
F sa , ^6:5 . F 6:2;3 . F i^,6 , which are depicted in fig. 8 and defined in eqs. (IV.10) and (IV.14). For 
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Ci 



the first two integral functions the pole comes from a tree propagator, while for the second two 
integral functions the loop integral itself generates the pole. (The F functions have the kinematic 
poles scaled out, but the original integrals defined in eq. (IV. 1) and related to the F in eq. (IV. 15) 
contain them.) 



7 ^6:2 
(a) 

6 1 



„lm 2 



(b) 



4-7 \ 
Q r 2m h 
^6:2;3 

(c) 



6 

A 1 



r 2m /i 
^6:2;6 

(d) 



[3] 



Figure 8. The four box functions which may appear as coefficients of a t 2 pole in a six-point 

[31 

amplitude; the dashed lines represent the t 2 channel. 

For the N = 4 supermultiplet with one gluon, four Weyl fermions and three complex scalars 
(with two states each) all ultraviolet and collinear singularities cancel leaving only the soft ones, 

A£r 4 (l+,2+,3+,4-,5-,6-) =-c r < ee (l + ,2+,3 + ,4-,5-,6-)i^f^^) . 

singular £ — * 



[31 

The tree amplitude in this equation contains a pole in the t 2 channel given by 



(4.5) 



Af ee (1+ , 2+ , 3+ , 4" , 5" , 6" ) = A%°° (2+ , 3+ , 4" , K~ ) A^ ee (K+ , 5" , 6" , 1+ ) + non-pole . (4.6) 

The four box functions which appear as coefficients of the P°l e i n the amplitude Aq.^ a have 
infrared singular behavior, 



7711m 
-^6:2 


singular 


r-ilm 
^6:5 


singular 


T?2mh 
r §:2;3 


singular 


T?2mh 
r 6:2;6 


singular 



l r 



l 



(-.s,,i) ^ f (-.s,o) ' - (-4 3] ) _e 



(-S 34 )- e + (-S 23 )- e " (~t 2 



l r 

2^2 

1 

2? 



(4.7) 



(-512)- 6 + 2(-4 3] )" e - (- S56 )- e " (-S34)- 

"(- S45 )- e + 2(-4 3] )" e - (-^ 23 )- e - (-aei) - 

and we have redundant constraints for determining their coefficients from eq. (4.5). This fixes the 
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relative coefficients of the four box functions to be unity so that 



^r 4 (l + ,2 + ,3+,4-,5-,6-) = c r ^ ree (2+,3+,4-,^-)- M < ee (K + ,5-,6-,l+) 



2 + + + F^) + non-pole , 



(4.8) 



where K = fc 2 + k 3 + k± and the tree amplitudes with the other helicity configuration of the 

[31 

intermediate line vanish. Thus, in the full amplitude, the t 2 pole must multiply the four box 
functions in fig. 8 with a relative coefficient of unity. 

From Table 2 we may also read off the factorization function in the channel to be 



^= 4 (4 3] ;fci,---,*6) 



-Sl2) ' - (-«56) £ - (-S34) 



1 (~g56) £ (-S34) 6 1 , 2 [ £12 

2e 2 (-*i 2 )" e + 2 m LP! 



1 ,_2 / s 34 \ 1 ,2 / «56 \ 7T 



- - In 



[3] 



+ {/c 2 <-> fc 5 , fc 3 <-> fc 6 , fc 4 <-> fei} . 



(4.9) 

f3l 

Note that the integral functions do not undergo any particularly large simplification as t 2 — ► 0, in 
contrast to the simplification appearing in the splitting functions. 

Finally, after explicitly taking the limit t 2 —> we may rewrite eq. (4.8) as 



1 



^r 4 (l + ,2 + ,3 + ,4-,5-,6-)^c r /. 2e < ee (2+,3 + ,4-,^-)^< ee (K + ,5-,6-,l 



2(F 0m (fc 2 , k 3 ,k 4 ,K) + F 0m (K, k 5 ,k 6 ,k 1 ) + = 4 ) 



(4.10) 



where the F 0m zero-external-mass box function is given in eq. (IV.lOf). The three terms in eq. (4.10) 
correspond to the ones in eq. (4.1), after re-expressing the functions in terms of the four-point [35,10] 
loop amplitudes. 

f3l 

We may obtain the pole in the t 3 channel by symmetry under reflection of the external states, 
so that the same results hold but with relabeled indices 

<f 4 (l + , 2+ , 3+ , 4" , 5" , 6" ) = c^A'r^A'r x 2 (f^ + + F*%* + F^) + non-pole . 

(4.11) 

Now compare the results in eqs. (4.8) and (4.11) against the explicitly computed result of 
ref. [11] 



<f 4 (i + ,2 + ,3+,4-,5-,6-) = Cf 



B x + B 2 W 6 (2) + B 3 



(4.12) 
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where 



B 1 



, ([12] (2 4) + [13] (34)) ([31] (16) + [3 2] (2 6)) {iff 
(1 2) (23) [45] [56] (tftf - s 12 s 45 ) (t [ ?kf ] - s 23 s 56 ) ' 



B 3 



B 2 




■j-i , 



■j'+i ' 



(4.13) 



and 



-ilm 
6:i+3 



-<2m h 
6:2;i+l 



+ F f 



-i2m h 
6:2;i+4 ■ 



(4.14) 



The B\ means complex conjugating spinor products in Bi, (kj) <-> [j k] without complex conju- 
gating factors of i, and the subscript j — ► j + 1 means applying a cyclic permutation of the six 
momenta {1, 2, 3, 4, 5, 6} — > {2, 3, 4, 5, 6, 1}. The Bi satisfy the condition 



The amplitude (4.12) clearly satisfies the multi-particle pole conditions described above: it does 

[31 [31 [31 

not contain a pole in the t[ 1 channel and the coefficients of the t 2 and t 3 poles match the ones 
in eqs. (4.8) and (4.11). Thus the amplitude (4.13) is consistent with the constraints imposed by 
multi-particle factorization, providing a further stringent check on its correctness. 

The same type of analysis may be performed for larger numbers of external legs. However, the 
singularities in e will not uniquely specify the coefficients of the boxes since the number of boxes 
proliferates. Nevertheless, all functions entering into the factorization function T n are uniquely 
specified from the singularities in e. 



In this section we prove the factorization properties of amplitudes described in the two previous 
sections. Our proof is based on identifying all potential poles in a generic amplitude that arise 
either from tree propagators or within loop momentum integrals. We use the integral reduction 
procedure [36,37,38,39,40,23], reviewed (and slightly modified) in appendix I, to make the poles 
from loop integrals explicit. One might expect that all non-factorizing contributions should be 
related to the infrared divergences; our method of proof makes this connection explicit. 

To summarize appendix I, any one-loop Feynman diagram may be expressed as a linear com- 
bination of integral functions with four or fewer legs multiplied by rational 'reduction coefficients'. 
In the standard integral reduction procedure one obtains coefficients of integral functions which 
can contain dependence on e. Since we are linking all non-factorizing contributions to infrared 
singularities, any e-dependence would have to be tracked because it could induce finite shifts in 
the splitting and factorization functions when multiplying divergent integrals. Rather than dealing 
with e-dependent reduction coefficients, we find it simpler to keep all e-dependence in the basis 



B 1 + B 2 + B 3 = 2A* ree (l+, 2+ , 3+ , 4" , 5" , 6" ) . 



(4.15) 



5. Proof of Universal Factorization 
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of integral functions. Thus, in the appendix we modify the usual reduction procedure in order to 
eliminate e-dependence in the reduction coefficients; through O(e ) this leads to the basis of integral 
functions: (a) D = 4 — 2e scalar bubbles, triangles and boxes, (b) D = 6 — 2e scalar bubbles and 
triangles, and (c) D = 8 — 2e scalar boxes. In the dimensional reduction or FDH schemes, where 
the number of states is always fixed at their four-dimensional value, there are no other sources 
of e-dependence. Thus, all amplitudes may be expressed as linear combinations of these integral 
functions with rational coefficients containing no e-dependence. 

If the reduction coefficients contain poles or if the integral functions contain discontinuities 
or poles, the amplitudes may not factorize in a simple way. After having performed the integral 
reduction procedure on all Feynman diagrams contributing to an amplitude, poles in kinematic 
variables arise from three sources:''' 

1) diagrams which have a pole coming from an explicit tree propagator, 

2) poles in the integral reduction coefficients, and 

3) explicit poles appearing within the box and triangle integral functions in the basis. 

The naively factorizing contributions are of the first type but all three can lead to non-factorization. 
We now systematically collect the non-factorizing pieces to the loop splitting function associated 
with each of the three types of kinematic poles. Subsequently, we will show that they are uniquely 
fixed by the singularities in e via eq. (2.6); this is because no non-zero linear combination can be 
constructed which is free of e _1 . 

5.1 Kinematic Poles from Tree Propagators 

Consider the Feynman diagrams containing a tree propagator which has a kinematic pole in 
the variable as depicted in fig. 5. In the limit that the kinematic variable vanishes one expects 
the amplitude to factorize on this pole since the diagrams break into a product of lower-point tree 
and loop amplitudes. However, if additional infrared divergences develop in the loop integrals as the 
kinematic invariant vanishes there will be discontinuities in the integral functions. In the off-shell 
case (t^ 7^ 0) there may be, for example, terms such as ln(— t^) which correspond to poles in e 
in the on-shell case (ij r ' = 0); in this case the off-shell to on-shell transition is discontinuous, and 
there is a 'discontinuity function' modifying the tree pole. 

To find the non-factorizing contributions in the loop integrals of the diagrams in fig. 5 (including 
the case of collinear limits), we categorize all discontinuities in the integral functions (coming 
from the reduction) as the kinematic invariant K 2 vanishes. (We denote momenta which arc 
sums of external momenta ki by an upper case letter; the lower case ki we reserve for on-shell 
external momenta.) Out of all possible integral functions appearing after the reduction, discussed 
in appendix I, only those where one leg has momentum K have a discontinuous off-shell (K 2 / 0) 
to on-shell (K 2 = 0) transition. We now step through the possible discontinuities coming from 
these integrals, which are collected in appendix IV. 

t We are assuming that the spinor helicity reference momenta are chosen so as to not introduce any additional 
poles in the channel of interest. 
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Consider first discontinuity functions arising from the D = 4 — 2e bubble function. Since 
bubbles are a function of a single momentum, the only ones which can have a non-smooth limit 
as K 2 — > are those with momentum K flowing through them. Bubble integrals with on-shell 
(massless) legs vanish by a standard dimensional regularization prescription [31], so the off-shell 
bubble (IV. 2) is itself a discontinuity function. Thus we have the discontinuity function given by 

b{ K 2 ) = ±I 2 [1](K 2 ) = -^L_(-K 2 )-* , (5.1) 

where r r is the ratio of Gamma-functions in eq. (IV. 3). 

Consider now the discontinuity functions arising from the D = 4— 2e triangle integrals, depicted 
in fig. 9. Starting from the three-external-mass integral, depicted in fig. 9a and explicitly defined 
in eqs. (IV. 7a) and (IV.4a), and taking K$ — > one obtains 

T^(K 2 , Kl Kj) T$™(KlKi) - d 2 (K 2 ;K 2 , iff) + d 2 (K 2 ; K*, K 2 ) , (5.2) 

where Tg m is the two-external-mass triangle depicted in fig. 9b and explicitly given in eqs. (IV. 7b) 
and (IV. 4b). Equation (5.2) defines the discontinuity function d 2 to any finite order in e. Through 
0(e°), it is 

MKl.KlKl) S ^(-K!)-' - ^f:^' - *(l - f ) ■ (5-3) 
and satisfies the property 

d 2 (Kl-KlKl) = -d 2 (Ki;KlKf), (5.4) 
as may be checked using the dilogarithm identity 

Li 2 (l -x) + Li 2 (l - x- 1 ) = ~ ln 2 (x) . (5.5) 




K 2 



K, 



\ 





K x K x K x 

(a) (b) (c) 

Figure 9. The triangle functions generate two types of discontinuity function as one of the 
external legs becomes on-shell (or massless). 



Next, by taking JTf — ► the two-external-mass scalar triangle in fig. 9b becomes 



= T^(K 2 1 )-d 1 (K 2 2 ), 



(5.6) 
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where the discontinuity function is 



di(K 2 ) = -^r(—K 2 )~ e . (5.7) 



Finally, by taking K\ — > for the one-external-mass triangle depicted in fig. 9c we obtain 

T^(K 2 )=d 1 (K 2 ). (5.8) 

Note that the zero-external-mass triangle vanishes. 

Now consider discontinuities in the D = 4 — 2e box. These are readily obtainable from the 
general relationship between box and triangle functions [40] 



h[l] 



E^^W + (-l + 2 e )c /f= 6 - 2 11] 



i=l 



(5.9) 



where 

4 4 



j = l i=l i,J = l 

Here pj = X^=i Kj-> where the iTj are the external momenta of the box functions and po = 0. 



The six-dimensional box function /f~ 6-2e [l] in eq. (5.9) (expanded to a finite order around 
e = 0) has a smooth transition between external massive and massless kinematics since it does 
not contain infrared divergences and the loop integral converges uniformly. Thus from eq. (5.9), 
the only discontinuity functions that may appear in box functions are linear combinations of ones 
appearing for triangle functions (5.7) and (5.3) (and explicitly listed in Table 5 in section 6) and 
there is no need for a separate analysis. 

Finally, consider the higher-dimensional integrals. These higher dimension boxes, triangles and 
bubbles asise from the reduction of tensor integrals, following the procedure discussed in appendix I. 
(If one were keeping terms of 0(e) and beyond, then the basis is further enlarged to include higher- 
dimensional integrals with n > 5 legs.) In general, these integrals are smooth as an external leg 
makes an off- to on-shell transition since no infrared divergences develop. Therefore they will 
not contribute a discontinuity function. There are, however, two exceptions which need a closer 
examination: 

[l{ ' 2e(l-2e)(3-2e) 1 1 

r lm,D=6-2er 11 / j^2\ _ r T I Tf2\-e 

/s 26(1 -6)(l-26) { ~ K) ■ 

These integral functions are proportional to I2 in eq. (5.1), but with e-dependent coefficients. The 
possibility of the integrals (5.11) contributing to non-factorization would be problematic since their 
difference is simply a constant, which would not be linked to a singularity in e. In appendix III 
we show that the integrals in eq. (5.11) do not contribute for K 2 — > 0, leaving only D = 4 — 2e 
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integrals as source of discontinuity functions; there are no additional discontinuity functions from 
any of the higher-dimensional integrals. 

Converting to the n-point kinematic variables in the K 2 = tf^ — > channel, the discontinuity 
functions which may enter from the loop diagrams in fig. 5 are thus 

d 2 it i ,t i+r ^ )__(-<.) - ^2 ( _ t [^'])- £ — L12 1 1_ ^fmJ' (5 - 12) 

°^ J e(l-2e) 1 tl j ' 

where 2<r'<n — r — 2. These discontinuity functions describe the additional contributions one 
would find from diagrams containing a tree propagator pole in For example, as K 2 = — ► 0, 
the first set of diagrams on the right-hand-side of fig. 5 behave as 

£ 4°+T(- • • >^ A ) ^2 (^" A > ■ ■ ■) + 4a E discontinuities , (5.13) 

A=± 

where the discontinuities are the ones in eq. (5.12) multiplying rational coefficients (determined in 
section 5.4) where the pole in K 2 comes from the intermediate tree propagator. That is, the dia- 
grams factorize as one might naively expect, except that there are possible additional discontinuity 
contributions multiplying the explicit tree pole. 

The origin of the discontinuities in eq. (5.12) may be traced back to the implicit expansion of 
amplitudes in e; if one were to keep the full functional behavior in e there would be no discontinuity 
for e < since they can be analytically continued to zero. Note that the discontinuity functions 
contain non-singular (in e) contributions, including the dilogarithm; these finite contributions are 
linked to the infrared singular ones because they are part of a single function. 

5.2 Kinematic Poles from Integral Reduction Coefficients 

As discussed in appendix II poles may appear in the reduction coefficients as two external 
momenta become collinear (or as an external momentum becomes soft) and not in multi-particle 
channels for general kinematics (except for the tensor bubble (1-20) which generates a multi-particle 
pole that is canceled by the D = 6 — 2e bubble). For convenience we choose the s\2 channel. 
Only loop diagrams of the type in fig. 10 may generate potential s\2 poles in the reduction 
coefficients. (For gluon amplitudes with scalar loops these types of diagrams were explicitly analyzed 
in refs. [13,15].) We will present a general argument and analysis that all such poles are proportional 
to the discontinuity functions (5.12). 




Figure 10. The class of loop integral which may generate a pole in S\2- 
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For the case of loop integrals which have a uniform non- vanishing mass in the loop (see eq. 1.1), 
it is a simple matter to argue directly from the loop momentum representation that there are no 
massless poles. The mass acts as an infrared cutoff which prevents the integral from diverging as 
a kinematic variable vanishes. Alternatively, in the integral reduction framework, the reduction 
coefficients explicitly given in appendix I do appear to contain the poles collected in appendix II; 
for the case of a uniform mass in the loop these poles are necessarily spurious (i.e., the residue 
vanishes) since massive loop integrals cannot contain massless poles. (The spurious nature of these 
poles was first noted by Brown [36].) The specific details of how the residues of poles in the 
reduction coefficients all vanish is in general quite intricate and may involve Abel's [41] dilogarithm 
identity; nevertheless they must vanish. 



As shown in appendix II, the reduction coefficients for the case of a uniform mass in the loop 
contain exactly the same apparent poles as for the massless case. In the massive loop case the 
residue of all massless poles cancel so one might expect the residues to also vanish for a massless 
loop. Indeed by taking the massive loop integral residues (which vanish) and taking m — > one 
obviously still obtains zero since the limit is smooth; all potential non-smoothness would be of the 
form ln(m), but such terms (as well as all others in the residue) vanish. 



We must, however, be a bit more careful for the massless case since we must take the limits in 
the reverse ordering: first we set m = and then we extract the residue of the kinematic pole. After 
setting the internal mass to zero, the extraction of the residue by taking the appropriate kinematic 
variable to vanish may no longer be smooth since new infrared divergences may develop. The new 
singularities can prevent the necessary identities for the vanishing of the residues in the massive 
case from holding. Thus we need to collect the new singularities and all associated discontinuities 
from the integral functions. This situation is completely identical to the case where the pole arises 
from a tree propagator. The analysis is completely identical and the non-smooth behavior of any 
integral function when a kinematic variable vanishes is described by the discontinuity functions 
(5.12). Thus any (non-spurious) poles in the reduction coefficients must give a contribution to 
splitting or factorization functions proportional to discontinuity functions. 



We have performed a number of explicit checks to verify the above general argument. Consider 
for example the scalar pentagon integrals with the kinematic configuration depicted in fig. 11, where 
momenta fei and are on-shell and the remaining Ki are either on-shell or off-shell. As discussed 
in appendix II, for this kinematic configuration the reduction coefficients contain poles in s\i- 
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^3 
K 5 

Figure 11. The original n-point integral may be reduced into a sum of pentagons; the type 
of pentagon depicted here is the only one which can generate poles in Si2- 

From the discussion in appendix I, the pentagon integrals are 

5 



4 = iE c ^4+ec / 5 D = 6 - 2e , (5.14) 



./W i r „_ T D=6-2e 
i=l 

where 

5 5 

Ci = J2 S ij 1 ' c = ^2a, Sij = -\{p i - 1 -p j - 1 f . (5.15) 

j=i i=i 

For on-shell k\ and k 2 it is straightforward to check that as s 12 — > 

^<»-<"' 

where 

4v / det5W = S23S34 - ssi^f , 4 Vdet S( 2 ) = s 34 s 45 - KlKl , 

4 V / det 5( 3 ) = 545^51 " S23^| , 4\/det S( 4 ) = S 5 1S12 , ( 5 - 17 ) 

4 V det 5( 5 ) = s 12 s 23 . 

The 5m„ are the matrices for the daughter box diagram obtained by removing the internal prop- 
agator between leg (i — 1) and leg i. (This follows the same labeling convention as for the box 
diagrams themselves.) Since Vdet 5< 4 > and Vdet 5( 5 ) are both proportional to s\ 2 , only the three 
coefficients ci,C2 and C3 contain a pole in S12 and are relevant. 

In all cases, the coefficients C4 and C5 do not contain any singularities. Rather, the integrals 
/| 4 ^ and themselves contain a divergence. Since we are only investigating the singular behavior 
of the reduction coefficients, we only consider the partial sum over ci, C2, and c 3 

^ + ^ + <^ ^ -T^ ^CTT g'-^ ^ - ^ + • 

(5.18) 

where we used eq. (IV. 8) to scale out the overall kinematic denominator from each integral function. 

As explained above, for any configuration of K 3 , K 4 and K 5 , being massive (Kf / 0) or 
massless (Kf = 0) we have 

lim (F 4 (1) - F 4 (2) + F 4 (3) ) = adi(si 2 ) + (3 d 2 (s 12 ; tff , s 45 ) + 7 d 2 ( ai2 ;tff , s 34 ) , (5.19) 
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where the a, (3 and 7 are constants. The combination of box functions and F^ is 

displayed in fig. 12. We have explicitly verified through O(e ) that this holds for all possible 
kinematic configurations. For each situation the real numbers a, (3, and 7 were found by explicitly 
performing the limit for all types of box integrals (i.e. for the different cases of whether or not 
K\ = 0, K\ = 0, or = 0). In all cases the discontinuity functions come from F^ 2 \ which is 
the only one of the three box functions which acquires new infrared singularities as S12 — ► 0. (See 
table 5 in section 6 for a list of all box discontinuities.) The identity (5.19) is quite non-trivial when 
actually performing the limit on the integrals in appendix IV since many dilogarithm identities must 
be used in the process. Nevertheless, in a theory with a mass such identities necessarily exist to 
cancel the pole in the reduction coefficients; in the massless limit these identities break down in the 
infrared, leaving behind discontinuity functions containing singularities in e. 
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Figure 12. An example of an non-trivial identity in the collinear limit k\ || k 2 is that this 
sum of box functions reduces to discontinuity functions; k\ and k 2 are on-shell while K3, K4 and 
are either on- or off-shell. 



As a second check, the general five-point integral with one power of loop momenta is [40,23] 



i=i j = l 



(j) fci||fe 2 



y^p*«(_i)<+i_ 

^ s 12 VdetSm VdetSW 



1YF, 



(5.20) 



so once again identity (5.19) guarantees that any pole in si 2 is proportional to a linear combination 
of discontinuity functions. One can continue in this way to show that any poles in S12 are propor- 
tional to linear combinations of discontinuity functions even for higher powers of loop momentum in 
the numerator. We have also verified that the poles in the reduction coefficients of scalar hexagons 
are also proportional to discontinuity functions. 

The explicit reduction results above are complete for an iV = 4 supersymmetric theory since 
the only integrals encountered [10] that can produce a collinear pole in reduction coefficients (using 
the reduction basis discussed in appendix II) are the ones discussed above. A proof of the lack 
of collinear poles in integrals for the case of external gluons with scalar and fermion loops has 
previously been constructed by a direct diagrammatic analysis [13,15]. Since gluon loops may be 
interpreted in terms of a linear combination of scalar, fermion and = 4 loops [1,42], this provides 
a complete check of our general argument for n-gluon amplitudes. 

In summary, all contributions with a kinematic pole in the reduction coefficients are propor- 
tional to 

y~] discontinuities , (5-21) 
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where the discontinuities are again the ones in eq. (5.12). The general argument we have presented 
above holds for any type of massless pole. Besides collinear poles, the reduction coefficients may 
contain poles as momenta become soft or ones not found in amplitudes, such as in non-adjacent 
kinematic variables (e.g. S13); in view of the above general argument all such poles must either be 
spurious or proportional to discontinuity functions. 

5.3 Kinematic Poles from Integral Functions 

The third source of kinematic poles are in the explicit forms of the integral functions contained 
in the reduction basis. In eqs. (IV. 7) and (IV. 15) of appendix IV overall dimensions have been 
explicitly pulled out from the D = 4— 2e scalar integral functions; the only ones that have kinematic 
poles in two- or multi-particle channel are I^ m , l\ m and J| mfl . In a two-particle channel s^j+i the 
integral functions 
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-.lm 



i+l 
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contain poles, which are depicted in fig. 13. These integrals arise from the reduction of the type of 
diagrams in fig. 10. All other scalar integrals in the basis do not diverge as s^j+i — > 0. We also note 
that the triangle function T^™ is equal to the discontinuity function di(si^ + i) defined in eq. (5.7). 
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Figure 13. The integral functions which have a pole in s^j+i. 
The case of multi-particle poles is similar, and only the two integrals (depicted in fig. 7) 
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contain poles in the channel t| r ' (r > 2). 
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Finally there are the higher-dimension integrals. Such integrals, however, cannot contain 
(non-spurious) poles. In general the nearest neighbor tf^ poles in the integrals come from isolated 
regions in the loop momentum integral. In six- or higher-dimensions the measure dPl prevents the 
occurrence of a pole from the dominant region of the loop momentum integral, where propagators 
become singular. This type of analysis is similar to the previously performed analysis of the 
suppression of poles in n-gluon diagrams containing internal scalars [13,15]. 

5.4 Linking all contributions to singularities 

Having found all possible non-factorizing forms containing poles in the channel of interest, 
we now go on to fix their coefficients by requiring that the singular terms are consistent with the 
known behavior in eq. (2.6). 

First consider two-particle collinear limits. The loop diagrams in fig. 2, containing the fac- 
torizing contributions to the loop splitting functions, may contain singularities in e which depend 
on the details (and gauge choices) of the process under discussion. A procedure which simplifies 
our discussion is to collect all singularities in the non-factorizing set. We therefore always subtract 
and add the discontinuity functions b(si^i) and c2i(sj,j+i) given in eq. (5.12) with coefficients ad- 
justed to move the singularities to the non-factorizing contributions. (See section 3.2 for an explicit 
example of this procedure.) 

In Table 3 we collect all non-factorizing contributions to Split loop in eq. (3.1), along with the 
singular terms which we use to fix the coefficients. These singular terms satisfy the property that 
there is no non-zero linear combination which is free of poles in e. Since all singular terms appearing 
in one- loop amplitudes are necessarily proportional to tree amplitudes, the final non-factorizing 
part of the loop splitting functions must be proportional to tree-level splitting functions, as given 
in eq. (3.9). 

To fix the coefficients we systematically step through Table 3. Starting with the first row we 
have the potential contribution to Split loop , 

^ ~(-zs iti+i y e -Li 2 (l-*) + 0(e), (5.24) 

which is the only non-factorizing function containing a ln(z)/e. We may fix the coefficient of this 
term by comparing it to the coefficient of the same term in the original n-point amplitude, on the 
left-hand-side of eq. (3.1). In the collinear limit where h = zK and fei+i = (1 — z)K, from eq. (2.6), 
we have 

-c r <sfcU a?" -crsS \y 2 +-M-ZT-) - \ ln ^l £ s P iit tree <T 1+ o( 

(5.25) 

where = (h-i + K) 2 . This is the only singular term in the original n-point amplitude 

containing a ln(z)/e in the collinear limit. 
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Table 3: Particular singular terms in the second column used to fix the coefficients of potential 
non- factorizing contributions to Split loop . 

Matching the coefficients of the ln(z)/e in the potential contribution (5.24) to the one in the 
n-point amplitude (5.25) gives the splitting function contribution 



c r 5ftSplit tree (d 1 ( Si , i+1 ) +^F^ +2 ) , (5.26) 

\ ki\\k i+1 / 



which is thus fixed by the coefficient S\_ 1 in the original n-point amplitude. Since F^f, 2 a ^ so 
contains a ln(— s i)i+1 )/e it is convenient to subtract out this this singularity using the discontinuity 
function (^(s^+i), so that we only adjust the coefficient of one singularity at a time. In the collinear 
limit the singular term (5.25) in the amplitude also contains a ln(— Si_i 5 x)/e which exactly matches 
a corresponding singularity in the loop amplitude ^4^°i on the right-hand-side of eq. (3.1). Pulling 
out the pre-factor c r Split tree , as in eq. (3.9), gives the first term in eq. (3.10). The collinear limit 
ki || ki + \ of eq. (5.26) gives the first entry of Table 1. 

Following similar reasoning, the coefficient of ln(l — z)/e in the amplitude matches the one in 
the splitting function 



Cr^SphtWdiK^) + ^F^ +3 ) , (5.27) 



where 



F% +3 ^ -1 (-(1 - z)s t , l+1 ) " - Li 2 (z) + 0(e) . (5.28) 

The contribution (5.27) corresponds to the the second entry in the third column of Table 2. 

The third and fourth rows in Table 3 contain singular terms which are nowhere to be found 
in n-point amplitudes or in the factorizing portion, so the coefficients must vanish, explaining the 
absence of such contributions in Table 1. 
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Since we have collected all ln(— s i)i+1 )/e singularities together by moving such terms from the 
factorizing diagrams to the non-factorizing contributions, the coefficient of the di^j+i) contribu- 

[nl 

tion is simply <5j to match the same singularity in the amplitude, as given in the third row of 
Table 1. Finally, the coefficient of ^(Si^+i), which is the only contribution to the splitting function 
containing a 1/e singularity (without a logarithm), is determined by the difference of the coefficients 
for the n-point amplitude and C[ n_1 l for the (n — l)-point amplitude. 

Now consider factorization in a multi-particle channel t\ — > 0. The procedure for adjusting 
the coefficients of potential contributions to the factorization functions is similar to the two-particle 
collinear case. Once again the non-factorizing contributions must be proportional to a product of 
factorized trees, because the singularities in e are necessarily of this form. Following the same logic 
the reader may step through the rows of Table 4 to obtain the entries in Table 2. 



Potential Contribution 


Selected Singular Part 


n:r— 


^■ln(-Si_ M ) 


Tp2m h 
n:n— r— l;i+r+l 


- ln(— Si +r -i^ +r ) 


J (Ar],[r'] Jr+r']\ 


- iH-t [ ;lit [ r' 1 ) 








l 

e 



Table 4: Selected singular terms associated with integrals with potential non-factorizing contri- 

\r] 

bution in the multi-particle t\ channel. 

Thus we have fixed the coefficients of all non-factorizing contributions in eqs. (3.1) and (4.1) 
for Split loop and J- n and established the rules in Tables 1 and 2. The discussion we have presented 
here is valid to any order in e; the only modification being that the box functions appearing in 
these tables be kept to higher order. 

6. Discontinuity Functions as a Tool for Evaluating Integrals 

In this section we illustrate the use of discontinuity functions as a tool for obtaining infrared 
divergent box integrals from known infrared finite ones. 

The infrared finite scalar box integrals of ref. [22] were evaluated in four dimensions. In 
ref. [23] a separate analysis for the infrared divergent integrals was performed due to the need for 
dimensionally regularized expressions. Since the discontinuity functions summarize the transition 
between infrared finite and infrared divergent integrals this suggests a procedure for obtaining 
infrared divergent scalar box integrals directly from infrared finite ones. After subtracting off the 
discontinuity function one simply takes the appropriate mass or kinematic invariants to vanish to 
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obtain the desired integrals. This procedure is made practical by eq. (5.9) which gives the necessary 
box discontinuity functions from the much simpler triangle functions. Thus, using known infrared 
finite box integrals, any infrared divergent box integral is reduced to an evaluation of simpler 
triangle integrals. 

To illustrate the method we reproduce all the infrared divergent box integrals (depicted in 
fig. 14) necessary for computations in massless theories given in ref. [23]. Starting with the four- 
mass box one takes the limit of a vanishing external kinematic invariant (or mass) and subtracts 
the appropriate discontinuity function. For example, to obtain the three-external-mass integral 
function (fig. 14b) from the four-external-mass integral (fig. 14a) given in ref. [22] we have from the 
first row of Table 5, 



lim 



F 4m {Ki ,K 2 ,K 3 ,K 4 )-d 2 {K( ■ Ki , {K x +K 2 ) 2 )-d 2 (K{ ; K{ , (K 4 +^i ) 2 ) 



F^(h, K 2 ,K 3 ,K 4 ) 
(6.1) 

This result agrees with the explicitly calculated result in ref. [23]. The reader may verify that all 
other infrared divergent box functions derived in ref. [23] may be obtained in this way, by using 
the discontinuity functions collected in table 5.* 

This method for obtaining box integrals may be used for any of the infrared divergent one-loop 
box integrals. In particular, box integrals with mixed massive and massless internal and/or external 
legs may also be obtained efficiently in this way. This method might also be useful for higher-loop 
integrals, but one would need a method for obtaining discontinuity functions that were obtained 
at one loop through eq. (5.9). 
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Figure 14. The box functions appearing in Table 5; the upper-case Ki represent off-shell 
(Kf ^ 0) and lower-case ki represent on-shell momenta. 



* In ref. [23] sign errors were inadvertently introduced in the four-external-mass box integral obtained from 
ref. [22]; these signs are corrected in ref. [11]. 
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Table 5: The discontinuity functions for all scalar box integrals encountered in computations in 
massless gauge theories. (See fig. 14.) 



7. Conclusions 

In this paper we have provided a general proof and discussion of factorization in massless 
amplitudes. Infrared divergences associated with massless theories significantly complicate matters 
as they induce contributions which are not interpreted directly in terms of a naive sum of products 
of amplitudes with a tree propagator. The loop integrals themselves may contain kinematic poles 
which cannot be interpreted in any simple way in terms of tree propagators. Furthermore, if 
the amplitudes contain logarithms in the kinematic variable in which the factorization is being 
performed, new infrared divergences develop. 

Nevertheless, we provided a proof that factorization of massless one- loop gauge theory ampli- 
tudes are described by a set of universal functions linked to the known [19,20,21] infrared diver- 
gences. In massless (or high-energy) QCD these singularities have been tabulated in previous papers 
[19,20,21]. The proof presented here was based on general properties of gauge theory amplitudes 
and not on a specific diagrammatic analysis. In particular, we made use of the reduction of any one- 
loop amplitude in terms of a basis of scalar box, triangle and bubble functions [36,37,39,38,40,23]. 
By identifying all possible sources of poles from the loop integrals, a proof was presented which 
links the coefficients of all such kinematic poles to the known infrared divergences appearing in the 
amplitudes. The non-factorization of amplitudes is described by a limited set of discontinuity and 
box integral functions, which were given explicitly. The non-trivial aspect of our proof is that all 
finite (in e) contributions to non-factorization are fixed in this way. 

In particular, we proved the conjecture in ref. [10] that the collinear splitting functions (or am- 
plitudes) determined from one-loop four- [5] and five-point [1,3] amplitudes are universal functions 
for an arbitrary number of external legs. We also presented a procedure for directly obtaining the 
splitting functions by calculating only two- and three-point diagrams. Although these diagrams are 
gauge dependent, after combining these diagrams with discontinuity functions whose coefficients 



33 



are fixed by the infrared divergences, the gauge invariant splitting functions are obtained. An ex- 
plicit sample calculation of the g — > gg splitting function in two different gauges was provided. We 
also calculated the contributions of massive fermions and scalars to these loop splitting functions. 

We have also shown how amplitudes factorize on multi-particle poles. Although there are 
contributions which cannot be interpreted in terms of simple factorization, these are given by a small 
number of functions whose coefficients are again fixed by infrared divergences. For calculations of 
six- and higher-point amplitudes multi-particle factorization provides a powerful constraint on the 
form of the amplitude. An explicit six-point example of the constraints imposed by multi-particle 
factorization was given. 

We have also outlined a procedure for obtaining the splitting functions to higher order in e. The 
higher order terms would be useful for performing phases space integrals at next-to-next-to-leading 
order using the formalism of ref . [20] . 

Although we have not explicitly discussed factorization as the momenta of external legs become 
soft, the same type of analysis performed here may be extended to that case. By a general argument 
presented in section 5.2, any kinematic poles in the reduction coefficients are proportional to the 
discontinuity functions. 

A spin-off from our analysis is an efficient method for generating dimensionally regularized 
infrared divergent box integral functions from the known finite massive box integrals. Indeed all 
box integral functions contained in refs. [23] may be obtained by taking the appropriate massless 
limits of the massive integrals after subtracting off an appropriate set of discontinuity functions 
determined by a recursion relation between between box and triangle functions. To extend this 
result to two-loops one would need a way of determining the discontinuities. 

The results presented here prove that universal factorization properties hold for any gauge 
theory one-loop amplitude with an arbitrary number of external legs. We expect factorization to 
continue to be a powerful tool in the calculation of gauge theory amplitudes. 

We thank L. Dixon, D.C. Dunbar, D.A. Kosower and A.G. Morgan for helpful discussions and 
suggestions. Research supported in part by the US Department of Energy under grant DE-FG03- 
91ER40662, in part by the National Science Foundation under grant PHY 9218990, and in part by 
the Alfred P. Sloan Foundation under grant BR-3222. 

Appendix I. Reduction of One-Loop Integrals 

In this appendix we review and slightly modify the reduction procedure [36,37,39,38,40,23] 
that rewrites n-point integrals in terms of lower-point ones for use in constructing our proof in 
section 5. This procedure allows any Feynman diagram to be expressed as a linear combination of 
m < 4-point integrals. The reduction procedure for tensor integrals is a bit different than for scalar 
integrals so we discuss these cases separately. 

1.1 Tensor Integrals 

First we review the reduction of Feynman integrals with a loop momentum dependent numer- 
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ator and with five or more external legs. Following this, the analogous decomposition of tensor 
box, triangle, and bubbles is presented. 

Consider then the reduction of any tensor integral function with more than four external legs 
(n > 4). Denote by I n [l ai ■ ■ ■ l aj ] the tensor n-point integral with j powers of loop momenta in the 
numerator of the integrand: 



/ B P«i...Z^]=i(-l)"+l(47T) 2 - e / 

J I'l 



A-2 H /Q1 . . . laj 



(2vr) 4 - 2e (Z 2 - m 2 )({l - Pl f - m 2 ) ■ ■ ■ {(I - p n _ x ) 2 - m 2 ) ' 

(LI) 

where the momentum routing is taken to be pi = Yl)=i-Kj, and the Kj are adjacent sums of 
the external momenta of the n-point amplitude. In gauge theory, in a Feynman-like gauge, the 
maximum number of powers of loop momentum in the numerator is n for an n-point function so 
j < n. Although we are interested in amplitudes with vanishing internal masses, we have inserted 
a uniform mass m in all internal propagators in order to display the difference in the corresponding 
reduction coefficients. We will show that the pole structure is independent of the mass, a fact we 
use in section 5 to explain why the poles in the massless reduction coefficients are proportional to 
discontinuity functions. 

The integrals in eq. (I.l) may be reduced by projecting the first component ot\ of the tensor onto 
a basis of four independent momenta; a simple choice of basis are the four vectors p\, P2, P3, Pa- 
The final results for the reduction of the amplitudes are independent of any particular choice of 
basis. Projecting the loop integral onto the four vectors yields 

I n [l a 1 P 2 . . . l a > } = pT A 1 + p^ A 2 +p^A 3 + p^ 1 A A , (1.2) 

where we have suppressed the indices a,i . . . ctj on the right-hand side of eq. (1.2) in Aj. The 
functions Ai are found by first contracting eq. (1.2) with the momenta pi, generating the four 
linearly independent equations, 



2I n [l- Pi l a *...n] = Y,UkA k , (z = l,2,3,4), (1.3) 

fc=i 

where 

t ik = 2pi • p k . (1.4) 

Since the products I ■ pi within the n-point integral on the left-hand side of eq. (1.3) may be 
expressed as 

2l- Pi = - [(I - Pi ) 2 - m 2 } + [I 2 - m 2 } + p 2 , (1.5) 

one can eliminate one of the Feynman denominators with each of the first two terms in eq. (1.5), 
which gives the integrand for an (n — l)-point integral. The third term in eq. (1.5) is independent 
of the loop momenta and gives an n-point integral with only j — 1 powers of loop momentum in 
the integrand. 



35 



Thus by solving the linear equations (1.3) for the Ai, we may express the n-point integral with 
j powers of loop momentum as a linear combination of n- and (n — l)-point integrals with j — 1 
powers of loop momentum. By inverting t we have 

Ai = j-^ In[2 i. Pk /-...r,] 

k T (1-6) 

= £ %^ (4- + i 1} r 2 ^ 3 • • • - • • • + vi^T^ ■ ■ ■ i a *i) , 

k=l 

where the elements [t]ki in (1.6) are the cofactors of t. The denominator A is the Gram determinant 
of the basis of four vectors, 

A = det(ttf) . (1.7) 

(Note the relative signs, between I n and I n -i, in the definition of the integral functions (1.1).) 

For dimensionally regularized expressions there is a technicality which must be addressed. 
Since the expansion on the right-hand side of eq. (1.6) is in terms of a basis of four-dimensional 
external momenta while the loop momentum on the left-hand side are (4 — 2e)-dimensional, one 
might worry that this could lead to an error. Consider, for example, an n-point integral with two 
powers of loop momentum in the numerator. The full (4 — 2e) expansion is 

I n [rf\ = £ p?p?A i3 + 5«%A , (1.8) 

where $"_2e] 1S a me t r ic which is non-zero only in the [— 2e]-dimensions. By contracting <5^ 2 e] 
eq. (1.8) we obtain 

I n [-l 2 } = -2eA , (1.9) 

where —l^ is the [— 2e] component of I 2 . After Feynman parametrizing the loop integral (1.9) and 
following the discussion in refs. [9], we break up the momentum integral as 

_ f d*i 4 d- 2 n t _ _ e i^y_ rsu (,,,-!-.. (L10) 



(2vr) 4 - 2e J (2vr) 4 - 2e T(l - e) J (2vr) 4 

Performing the integrals one finds 

I n [l 2 } = -el^= e -^[l}, (1.11) 

where I® =6 ~ 2e [l] is the n-point scalar integral in D = 6 — 2e dimensions (given in appendix IV. 4). 
Note that I^ =6 ~ 2e is completely finite for n > 3, since there are no infrared or ultraviolet divergences 
in six dimensions for n > 3. Thus 5^ 2< ^A in eq. (1.8) is of 0(e) and therefore does not enter into 
the reduction procedure through O(e ). One can easily generalize this argument to show that the 
difference between a reduction in a four-dimensional basis and a (4 — 2e)-dimensional basis is of 
O(e) as long as no ultra-violet divergences occur in the integrals. In renormalizable gauge theories 



36 



(for a Feynman-like gauge) all diagrams with five or more legs are ultra-violet finite, so that the 
reduction coefficients for n > 5 do not depend on e. (If one were to keep terms of 0(e) or higher, 
then one would need to keep the higher dimension integrals for n > 5, but these do not possess 
infrared divergences, discontinuities as an external mass vanishes, or massless poles. Following the 
discussion in section 5, these do not contribute to non- factorization.) 

For n < 4 the reduction procedure is similar except that the integral does not contain a 
complete set of external momenta that spans four dimensions, so we must include the metric 
tensor. The metric tensor must be taken to be 4 — 2e dimensional since ultra-violet divergences are 
encountered: boxes with four powers, triangles with two or more powers of momentum, and bubble 
integrals are all ultraviolet divergent. The reduction coefficients of such integrals can be expected to 
contain e-dependence. (In the case of supersymmetric theories, one can arrange all n > 3 integrals 
to have at least two powers less of loop momentum in the numerator so that all integrals with 
n > 3 are ultra-violet finite; this means that the reduction coefficients are independent of e in 
supersymmetric theories [11].) 

Consider the standard expansion of the tensor integrals for n < 4, 

I n [l a ] = A 1;iP f , (1. 12a) 

I n [l^n = A 2 , S^l] + ^pTp? , (I-12b) 

4[rrr] = a 3;i te e] P? 8 + ^XfT + <ff 4 -VH + a ^pTp?p? , (i-i2c) 



i n [i ai i a n a3 i ai ] =a 4;0 

+ A^-ij 



°[4-2e]°[4-2e] ~T °[4-2e] °[4-2e] ^ °[4-2e] °[4-2e] 



^X\pTpT + g y m - + ^kipTpTpTpT , (i-i2d) 



where the repeated indices are implicitly summed over and some of the coefficients are related by 
symmetry. By dotting the momenta pi into these equations and contracting various indices one 
obtains a complete set of equations which may be solved for the coefficients A. The equations 
obtained by contracting indices contain explicit e-dependence through ^_ 2e ] ^ap 2 ^ = 4 — 2e. 

We now outline a modified procedure which eliminates explicit e-dependence from the equations 
obtained by contracting indices. Consider the first integral (1. 12a); in this case the metric tensor 
does not appear and the Ai ;i may be solved by dotting external momenta into eq. (1. 12a); no 
e-dependence enters in these equations. 

Next we have the integral (1. 12b), which has a ^^e] i n the tensor expansion. To solve for 
the A2 ; o we perform instead the contractions with the tensor which is non-zero only for the 

[— 2e] dimensions. Performing the contraction of eq. (1. 12b) with S^_^ we obtain 

I n [/ £ 2 ] = -e/f= 6 - 2e [l]=2e^ 2;0 , (1.13) 

so that 

^2;0 = -^„ D=6 - 2£ [l]. (L14) 
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The coefficient of I,f~ 6 2e [l] does not contain e-dependence. The remaining equations for A 2;i j 
obtained by dotting into external momenta are manifestly free of e. 

For I n [l a H a2 l a3 ] we may contract eq.(I.12c) with 25^p% to obtain 

n-l 

2I n [lll-p 3 }=2eY J A^t ]i . (1.15) 

i=l 

Using the expansion 2l-pj = — [(/ — pj) 2 — m 2 ] + [I — m 2 ] +p 2 and eq. (1. 11) we obtain the expanded 
form 

2I n [l 2 l. Pj ] = -c(e + 1 1) ' D=6 - 2£ - tilf =6 " 2e +^ D=6 - 2£ ) • (1-16) 
By inverting eq. (1.15) we then have 

1 n— 1 

Note again that there is no e-dependence in the coefficients of the integrals. Further equations 
which are obtained by only dotting external momenta into eq. (1. 12c) are manifestly independent 
of e. 

Finally, consider the case of four powers of loop momentum in the numerator. If we contract 
eq. (I.12d) with 5^*5^* we obtain 

I n [lf} = -6(1 - e)I° = 8 " 2e = (4e 2 - 4e)A 4;0 , (1.18) 

so that 

A,, Q = \lZ=*~ 2 * . (1.19) 

Once again all e-dependence cancels from the coefficients. The other coefficients may be obtained 
by continuing the reduction process and are also free of e. 

Thus, by introducing the D = 6 — 2e and D = 8 — 2e dimension scalar bubble, triangle and 
box functions (given in appendix IV.4) we may avoid e-dependence in the reduction coefficients; 
this dependence has been pushed into these integrals. (In the dimensional reduction [28] or FDH 
[5] schemes there are no other sources of e dependence, since the numbers of particle states are 
the four dimensional ones.) Upon rewriting the higher-dimensional integral functions in terms of 
D = 4 — 2e ones, the conventional reduction (through C(e )) is regained. 

As an explicit example of the modified reduction procedure consider the tensor bubble integral 
in a massless theory 

h^V] = -i{^) 2 ~ e J 



d 4 ~ 2e i m v 



(2vr)4-2, P{l-K) 2 (1.20) 
= B 1 5^_ 2e] K 2 + B 2 K^K". 

First consider a conventional reduction procedure. By either tracing over the indices with 5^_ 2< ^ 
or dotting into the K^K V we obtain the two equations 

(4-2e)B 1 + B 2 = 0, {B 1 + B 2 ) = J/ 2 [l] , (1.21) 
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where we used 2K ■ I = —(I — K) 2 + I 2 + K 2 and also dropped the tadpole diagrams, which vanish 
in dimensional regularization. (For higher-point functions it is more convenient to dot only one 
momentum at a time into the expansion for the integrals, followed by iterating the equations for 
each power of loop momentum.) Solving the two equations we have 

B > = -l3^W- * = 5&N- <'- 22 > 

which contain explicit e-dependence in the coefficient of the scalar bubble integral; it is this e- 
dependence which we wish to avoid by using the modified procedure. 

Now consider the modified procedure, where we extend our basis of integral functions to include 
the six-dimensional bubble. We obtain the first equation by tracing over the [— 2e] dimensions in 
eq. (1.20) to yield 

I 2 [l 2 } = 2eB x K 2 

(1.23 

= -eI 2 D =^[l], 

which we may use to solve for B\. The second equation is obtained by dotting K^K U into eq. (1.20) 

Si + B 2 = i/ 2 [1], (1.24) 

so that 

* = -^= 6 " 2e > B2 = \m + ^^ 6 " 2e • 

Note that there is no e-dependence in the coefficients of the integral functions in eq. (1.25). It is 
easy to verify that this solution is identical to the one in eq. (1.22) after using eq. (IV.20) to express 
the D = 6 — 2e integral functions in terms of D = 4 — 2e ones. 

1.2 Scalar Integrals 

Consider now n-point scalar integrals with n > 5. Denote by I n \\\ the scalar n-gon integral 
with no powers of loop momenta in the numerator of the integrand: 

In[l] = <-l) n+1 (4vr) 2 e J (27r)4 _ 2£ (/2 _ m 2)(( Z _ pi )2 _ m 2) . . . ((/ _ ^2 _ m 2) ' ( L26 ) 

A convenient method for reducing the integral in eq. (1.26) is given in ref. [38]. This method is 
based on the observation that for four-dimensional momenta defined by pi = Yl]=\Kj, and for 
n > 5 point integrals a solution to 

6 6 

5>P? = 0, J> = 0, (1.27) 

i=l i=l 

can be found for some constants 6j. At least six 6, are required for a non-trivial solution since 
five or more vectors are linearly dependent in four-dimensions and the second equation provides an 
additional constraint. For n = 5 this solution breaks down because there are only four independent 
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momenta present within the integral; another technique [39,40,23] may be used to reduce scalar 
pentagons whose results we quote below. 

The reduction for n > 5 proceeds by first multiplying the integrand of eq. (1.26) by unity in 
the form, 

and using the properties of hi (i.e., Y2i^ 2 = and • Pi = 0) so that we may express the 

numerator in eq. (1.28) as 

6 6 

6i(p? -™ 2 ) = £ - K) 2 " ™ 2 ) , (1-29) 

i=l i=l 

which corresponds to a sum over factors in the denominator of the integrand in eq. (1.26). The mass 
dependence in the denominator of eq. (1.28) also drops out via Ei=i hm 2 = 0, and the original 
integral in eq. (1.26) may then be expanded as, 



m 2 ] 



^n[l]= 6 \ 2 Y,h I n [{l - Pi f - 

U ''''' <=J (1.30) 

=-^r 1 r ^E^e + 1 i) [i], 

where i^'li is the (n — l)-point integral obtained from the n-point integral by removing the prop- 
agator in eq. (1-26) between legs i — 1 and i. 

The solution of eqs. (1.27) are ratios of kinematic determinants, given in detail by Melrose [38] 
The result in solving for the bi is that the coefficients in front of the in eq. (1.30) are inversely 

proportional to 

det 



hp^-pj^f] . (1.31) 



2 

Any poles in the scalar integral reduction coefficients (for n > 6), for the case of a uniform or zero 
internal mass, must come from this determinant. Since the mass dependence m in the integral 
reduction coefficients drops out for the case of a uniform mass around the loop, the poles for m = 
orm/0 are the same. 

This leaves us with the special case of n = 5, that is, the reduction of scalar pentagon integrals 
down to boxes. For this case, we quote the results of ref. [40]. The scalar pentagon with a uniform 
internal mass and external momenta Ki may be expressed as 



^SQf+^f", (L32) 



i=l 

where 

5 5 

c i = '^2 S ij 1 > c = ^Q, = m 2 - \{pi-i -Pj-if ■ (1-33) 

j = l i=l 
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In summary, by iterating the above tensor and scalar reduction procedures from n- to (n — 1)- 
point integrals any D = 4 — 2e one-loop amplitude may be expressed as a linear combination of 
scalar bubble, triangle, and box integral functions. By using a slightly modified Passarino-Veltman 
reduction of tensor integrals we avoid all e-dependence in the reduction coefficients at the cost of 
introducing higher-dimensional scalar integrals. The explicit forms of scalar functions which may 
appear in massless gauge theories is summarized in appendix IV. 



Appendix II. Kinematic Poles in Reduction Coefficients 

In this appendix we find the potential kinematic poles in a particular channel within the 
integral reduction coefficients, reviewed in the previous appendix. We will systematically step 
through and discuss the various kinematic denominators of the coefficients to exhibit their structure. 
Furthermore, we find that one can side-step the poles in a given channel for the reduction of integrals 
down to n = 5 and n = 6 for tensor and scalar integrals, respectively. We also track the difference 
in the kinematic denominators when there is a uniform internal mass in the loop and when there is 
no mass, which is used in section 5 to understand why these poles are proportional to discontinuity 
functions. 

II. 1 Poles from n > 6 integral reductions 

The coefficients which may appear, in each step of the reduction of integrals with six or more 
external legs, are proportional to the inverses of the two kinematic determinants given in eqs. (1.7) 
and (1.31), 

(a) det[2p CTi • p a .] , (i, j = 1, . . . , 4) , 

1 (II.l) 

(b) det[--(p ai -p aj ) ] , = 1, ...,6), 

where the p <Jj (with j = 1 . . . n — 1) are any of the pi = J2j Kj- The Kj are external momenta of 
the integral function being reduced and may be on- or off-shell. The final result for the amplitudes 
is independent of any particular choice of the bases {p aj } at each step of the procedure. (This 
clearly holds for tensor reductions [36,39] and has been explicitly shown by Melrose [38] for the 
case of scalar reductions.) The first determinant (a) is the Gram determinant arising from tensor 
reductions, while the second (b) is a modified Cayley determinant coming from the reduction of 
scalar integrals. 

Whenever either of these determinants has zeros the reduction coefficients can have poles. The 
Gram determinant will vanish whenever the momentum basis of the reduction collapses, which 
may happen if two (even non-adjacent) massless external momenta are collinear or one external 
momentum is soft. Note that as a multi-particle kinematic variable vanishes, the two determinants 
considered as general analytic functions will not vanish (except for special kinematics, such as when 
all kinematic variables are time-like); thus we only need consider collinear (and soft) poles. 

In any channel where the momenta of given pair of external legs become collinear one can avoid 
zeros in the determinants (a) and (b) by a judicious choice of the momentum basis. For example, 
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by not including either p\ or p 2 in the basis, as the momenta of legs 1 and 2 become collinear the 
Gram determinant (a) does not vanish. (For n > 6 there are always at least five momenta from 
which to choose the four needed in the Gram determinant.) In this way we may avoid collinear 
poles in any given channel coming from the tensor reduction coefficients for n > 6. The same type 
of reasoning allows one to avoid zeros in the second determinant (b) for n > 6. For the case of 
n = 6 one may explicitly check that there are zeros. Thus by making use of arbitrariness in the 
momentum basis choice of the reduction one can avoid the poles in the coefficients in any channel 
of the n goes to n — 1 reduction, down to scalar hexagons (n = 6). The poles in the scalar hexagon 
reduction coefficient do not, however, depend on the uniform internal mass. 

In another channel one would choose a different momentum basis to again avoid poles in n > 6 
reduction coefficients. The consistency of this procedure follows from the fact that all final results 
for amplitudes are independent of the choice of basis in each step; a change of momentum basis 
only shifts the poles to different steps of the reduction. Our choice is to always push the collinear 
pole in any given channel to the scalar hexagon and the pentagon reduction, discussed below. 

II. 2 Poles from n < 5 integral reductions 

Now consider the potential kinematic poles that arise from the reduction of n < 5 point 
integrals. As discussed in the previous appendix, the only denominators which may appear are 



where m is a uniform mass in the loop. The first one comes from tensor reductions with n < 5 
and the second from reducing scalar pentagons as given in eq. (1.32). There are no multi-particle 
zeros in either determinant for n > 2. For n = 2 the multi-particle pole appearing in the reduction 
coefficient of the tensor bubble (1.25) multiplies I^ =6_2e , given in eq. (IV.20), which then cancels 
the pole. One cannot avoid collinear poles for n < 5 since there is only one basis choice for the 
reduction of these integrals; the poles which have been side-stepped above for n > 6 appear here. 
For example, as s 12 — ► the diagram depicted in fig. 11, with on-shell legs k\ and k%, is the only 
type of pentagon which has a pole in this channel. 

The Gram determinant (a), which has no dependence on the uniform internal mass, contains 
zeros whenever any pair of on-shell momenta Ki = ki become collinear. (Recall that lower case 
momenta denote on-shell external momenta of the amplitude.) 

For the case of m = 0, it is straightforward to investigate zeros of the second determinant (b) 
in multi-particle and collinear channels. As a multi-particle kinematic variable (r > 2) vanishes 
detS will not vanish (except for special kinematic configurations). There will however be zeros in 
detS for collinear massless momenta. For example, the kinematic configuration in fig. 11 gives 

det S = --^si2 (s23 S34 S45 S51 - «23 S12 S51 Kl + s 5 i K\ K\ S23 - S34 K\ S23 - K\ s 45 s\\ , 



(a) A = det[2pi ■ pj] = det [2^ ■ Kj] , 

(b) det S = det m 2 - -(pi_i - Vj-\f , 



(l<*,j<n-l<4), 



= 1, ... ,5) 



(II.2) 



(II.3) 
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which clearly vanishes in the limit s i2 — > 0. Thus, poles of the form may develop in the 

reduction of scalar pentagons to scalar boxes. 

For the case of m / 0, the second determinant (b) contains explicit dependence on the uniform 
internal mass. However, this dependence is relatively simple since 



detS 



= detS 

m/0 



+ —A . (II.4) 

m=0 16 ^ ' 



The Gram determinant has zeros whenever det S\ m= o has one, so that det S\ m ^o will have the same 
zeros as det S\ m= o- 

In summary, in any given multi-particle or collinear channel the integral reduction may be 
performed so that there are no poles in the coefficients for n > 6. For n < 6 poles may appear as 
momenta become collinear (or soft). However, one finds the same apparent poles in the reduction 
of loop integrals containing a uniform internal mass as for the massless case. 



Appendix III. Higher Dimension Integrals: Special Cases 



In this appendix we show that the two special cases, the higher dimension integrals I®~ 6_2e and 
/j 1 " 1 ' D — 6_2e [ n e q_ (5.n) ; do not contribute as discontinuity functions. Both are proportional (with e 
dependent coefficients) to the D = 4 — 2e scalar bubble in eq. (IV. 2), which becomes a discontinuity 
function as the single external mass vanishes. The potential appearance of these integrals as 
discontinuity functions might introduce contributions to the splitting and factorization functions 
not linked to the singularities in e. In this appendix we show that such potential contributions do 
not occur. 

For both the bubble /^ >=6_2e (K 2 ) and triangle / 3 lm ' D ~ 6 ~ 2e (K 2 ) we will examine the appear- 
ance of poles in the K 2 channel. Poles may in principle come from either tree propagators, as dis- 
cussed in section 5.1 or from reduction coefficients, as in section 5.2. (These two higher-dimensional 
integral functions themselves do not have poles.) 

First consider the appearance of 6_2e . As discussed in appendix I, the Passarino-Veltman 
reduction takes an n-point integral function with m powers of loop momentum and reduces it to 
a combination of n and n — 1 point integral functions with m — 1 powers of loop momentum. 
Thus the only way to obtain a tensor bubble (two-point) integral function with two powers of loop 
momentum, which generate the D = 6 — 2e scalar bubbles via equations (1.20) and (1.25), is to start 
with an m-point integral with m powers of loop momentum in the numerator. Thus, to obtain 
D = 6 — 2e bubble functions we must consider (maximal) tensor integral functions of the form 

4 * — — o , (ni.1) 

(2vr) 4 - 2e I 2 (I - i^) 2 • • • (/ - K m f 

In ordinary or background field Feynman gauge this maximum occurs for diagrams where all legs 
attached to the loop are gauge boson legs [11]; diagrams where fermion lines are entering or exiting 
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the loop will have one less power of loop momentum. Thus we need consider loops with only gluons 
attached. 

First consider the case where the kinematic pole arises from a reduction coefficient. As dis- 
cussed in appendix II, no multi-particle poles arise from reduction coefficients. For the case of 
two-particle (collinear) poles one can side-step the poles from tensor integrals down to pentagons, 
as discussed in appendix II. 1. Furthermore, for the pentagon and below, using the formulas of 
appendix I, it is straightforward to check that the coefficients of the D = 6 — 2e bubbles contain at 
most a single pole in any given channel. Since /| )=6_2e [l](i^ 2 ) is proportional to (K 2 ) 1 ~ e the pole 
is spurious since it is canceled. Alternatively, one may investigate the m-point loop integrals with 
m powers of loop momentum directly; this type of analysis has already been discussed in ref. [13,15] 
with the result that they do not contribute a pole to the amplitude. 

Now consider the bubble function J^ )=6_2e from diagrams with a tree propagator as on the 
right-hand-side of fig. 5. The bubble functions can arise either directly from diagrams with two-point 
loops or from the integral reduction of higher-point loop diagrams. The diagrams with a bubble 
loop are part of the 'factorizing' contribution and therefore do not concern us in this appendix. For 
the bubble functions that arise from higher-point diagrams we show below that all integrals with a 
maximum number of powers of loop momentum suppress the pole from the intermediate leg, since 
they are proportional to either 

K 2 5^X] or Rai > ( m - 2 ) 

where the index a\ dots into the gluon (or gauge boson) line of the intermediate factorized leg and 
Ki = K in eq. (III.l). The integrals cannot produce any further poles in external mass K 2 , as 
can be seen from the fact that such poles do not exist in the reduction coefficients, discussed in 
appendix II, or in the basis of integral functions, discussed in appendix IV. (The D = 4 — 2e single 
mass triangles are an exception, but such integrals with more than four-point kinematics do not 
occur in the reduction of diagrams with tree poles that we are considering.) The first type of term 
is clearly sub dominant in the factorization since the factor of K 2 cancels the pole from the tree 
propagator. The second type of term is also subdominant because the gluon contracts against a 
tree amplitude which makes up a conserved current; the longitudinal term vanishes sufficiently fast 
to be irrelevant. 

First we contract eq. (III.l) with the external momentum K^ 1 and rewrite the scalar product 
as 21 ■ K\ = — (/ — Ki) 2 + I 2 + K 2 . This produces three new integrals, 

2I m [l ■ K\ I " 2 . . . |"»] = /^[P 2 . . . l am ] - i£-i[* aa • • • l am ] + K 2 I m [l a2 . . . r m ] , (III.3) 

where I^-i is the loop integral with m— 1 legs obtained from (III.l) by removing the jth propagator. 

The first two integrals are functions of Ki + K 2 and K\ + K m and will never reduce further 
into a bubble integral with momentum K" flowing through. The third integral does not possess a 
K 2 pole within the integral function and is then suppressed by a factor of K 2 in front. Thus in the 
factorizing limit there is no unsuppressed contribution. 
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Dotting the integral (III.l) into Kf leads only to potential discontinuities suppressed by Kf; 
this necessarily means that the integral (III.l) must be proportional to either of the two tensors in 
eq. (III. 2), which as discussed above are suppressed in the factorization limit. 

The single-external-mass triangle I^™ 1 ' D ~ 6_2e (IV. 21) also requires special care as it is discon- 
tinuous in the limit that the external mass vanishes. The argument of these integrals are the s^j+i 
where i and i + 1 label the momenta of two massless legs of the triangle. Such triangles are only 
discontinuous for the collinear or soft limits of ki and fcj+i. In the integral reduction, this triangle 
function can appear from two places. The first is from the factorizing diagrams displayed in fig. 2. 
The integrals from these diagrams are directly taken into account in our analysis since we calculate 
them explicitly; indeed, a combination of D = 4 — 2e and D = 6 — 2e integrals from these diagrams 
is what makes Split fact non-zero in eq. (3.7). 

The second place / 3 lm ' D - 6 ~ 2e mav arise is in the reduction of integrals of the form in fig. 10, 
which is given for the S12 channel (all remaining Sj^+i channels are similar); in this figure Ig" 1 ' D - 6_2e 
would arise in the reduction term where all loop propagators between and k n (following the clock- 
wise ordering) are removed; the single external mass in this case is Si2- As discussed in appendix II, 
down to pentagon integrals we can avoid poles from tensor integrals in any given channel by choos- 
ing an appropriate momentum basis in which to perform the reduction. For the pentagon to box, 
and box to triangle reductions it is straightforward to verify, using the formulas of appendix I, that 
the coefficient of / 3 lm ' D ~ 6_2e (si2) does not have a pole as S12 — ► 0. Thus / 3 lm ' D - 6 " 2e (si 2 ) does not 
contribute to the non-factorizing parts of the splitting and factorization functions. 



Appendix IV. The Basis of Integral Functions 

In this appendix we collect the integral functions useful for the discussions in the text; these 
integral functions were obtained from ref. [40,23]. The n-point scalar one-loop integral in 4 — 2e 
dimensions is 

d^- 2e p 1 



I n = (-l) n+1 i{^f- i J 



( 2vr )4- 2e p 2( p _ K j2(p -K 1 -K 2 f---{p-K 1 -K 2 K n _ x , 



)2 > 
(IV.l) 

where K i: i = 1, . . . , n are the external momenta, which may be either on- or off-shell. As discussed 
in appendix I any gauge theory one-loop amplitude can be reduced to a linear combination of (a) 
D = 4 — 2e scalar box, triangle and bubble integrals, (b) D = 6 — 2e scalar triangle and bubble 
integrals and (c) D = 8 — 2e box integrals. The higher dimension scalar integrals are defined by 
replacing the 4 — 2e in eq. (IV.l) with the appropriate dimension D. Following the conventions of 
ref. [23] integrals without a dimension label are taken to be in D = 4 — 2e. We first present the 
D = 4 — 2e integrals and subsequently give the higher dimension ones. 

IV.l Bubble functions 

The D = 4 — 2e two-point integral function is 

« K2 » = iT^<- K2 »-' (IV - 2) 
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where 



r(i + e)r 2 (i- e ) 

r(l - 2e) 



(IV.3) 



IV. 2 Triangle functions 

In fig. 9 the three types of triangle integral functions that may appear in massless gauge theories 
are given. The three D = 4 — 2e scalar triangle functions are 



lt(Kf,Ki,Ki: 



— Y 



Li 2 



1 + iSj 
1 - iSj 



- Li 2 



1 ~ iSj 
1 + id. 



ir(Kf,Ki) = 



(-Kf) - {-Kl) 



-2\-l-e 



where 
Si 

and 



a/a; 



5o = 



-K\ + Kl - Kl 



A 3 = -{Kl) 2 - (Klf - {Klf + 2K{Ki + 2K%K% + 2K'£K't . 



-2\2 



-2\2 



2 T^2 



'2 7^2 



-^l 2 - gj + gg 

a/a; 



'2 7^2 



(IV.4a) 

(IV.4b) 
(IV.4c) 

(IV.5) 
(IV.6) 



For our purposes it is a bit more convenient to deal with integral functions where the denom- 
inators have been scaled out so we define 



l!-(Ki,KlKi: 



rr(Kf,Ki,Ki: 



li™(KlKi) = Y^TJqT^iKl Kl) , 



ll m {Kl) 



(IV. 7a) 
(IV 7b) 
(IV 7c) . 



IV.3 Box Functions 

The scalar box (four-point) functions appearing in computations with massless internal lines 
have already been extensively discussed in refs. [22,23]. The reader is referred to these papers for 
further details. It is convenient to define this function as 



F(K 1 ,K 2 ,K 3 ,K 4 ) = -^*lh, 



(IV.8) 



where the symmetric 4x4 matrix S has components (i, j are mod 4) 



Si3=--(Ki + --- + Kj-i) 2 , i^j; S u = 0. 



(IV.9) 
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The external momentum arguments K\ mmm ± in equation (IV.8) are sums of external momenta ki 
that are the arguments of the n-point amplitude. From ref. [23], through O(e ), we have (after 
correcting a sign in the first box function) 

F 4m (K 1 ,K 2 ,K 3 ,K 4: ) = i/-Li 2 (|(l-Ai+A 2 + p)) + Li 2 Q(l-Ai + A 2 -p)) 
" Li 2 (-5Sr(l " Ai - A 2 - p)) + Li 2 (-^(1 - Ax - A 2 + p)) 



2 \X 2 J \l + X 1 -X 2 -p 



(IV. 10a) 



F 3m (k u K 2 ,K 3 ,K 4 ) = -\ 



1 H-K 2 )(-K 2 ) 
2e 2 V H) 



+ (~t)- e ~ (-Kir- - (-Ki)-' - (-Kl 

i ((-KD(-Kiy 



2e 2 



(s) 



Lio 1 



K 2 K 2 



st 



F^(h,k 2 ,K 3 ,K 4 ) = -- (s)-* + (-t)^ - (-Ki)-^ - (-Ki) 



J )+^ a (f)« (IV - 10b) 
~Kl)(-Kl 



1 

2^2 



(-*) 



+ Li 2 1 



Kd\ ( K 2 \ 1 2 ( s 

T +Ll2 + 2 ln \* 



(IV.lOc) 



-2\-e 



F (k 1} K 2 , k 3 , if 4 ) = --[(-^)- £ + ("i)" e " (-#2)" e " (-^4) 
F lm (k 1 ,k 2 ,k 3 ,K A ) = -1 + (-t)" e - (-K|)- e 

+Li2| £ 1 _M) +Li2 ( 1 _^ + ^ ( £)^ 



^(^,^^3^4) = -^ (-s)" £ + (-t) 



1 , 2 71-2 

+ - In 2 



2 ' 



(IV.lOd) 

(IV.lOe) 
(IV.lOf) 



where the fej denote on-shell momenta and the if, off-shell momenta. The kinematic variables are 



s = (h+k 2 ) 2 , 



t = (k 2 + k 3 f 



(IV.ll) 



or with k relabeled as K for off-shell (massive) legs and the functions appearing in F£ m are 



p = Jl-2X 1 -2X 2 + Xj-2X 1 X 2 + X 



2 > 



(IV.12) 



and 



Kl Kl 



Ai = 



(Ki + K 2 ) 2 (if 2 + K 3 ) 2 ' 



A 2 = 



(Ki + K 2 ) 2 (K 2 + K 3 ) 2 



(IV.13) 
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When dealing with the kinematics of an n-point amplitude it is convenient to label the integral 
functions in terms of the kinematic variables which may appear. Thus, we define 



n4m 



n:r,r' ,r 



F {ki + • • • + ki+ r —i, kj,-\- r + • • • + hi-\- r -\- r i — i , + • • • + A^_|_ r _|_ r 



ki+r+r' +r" + " " " + l) j 



/_i_ r //_l, 

(IV. 14a) 



n2m h Z?2m h 



F {ki — 2i f^i — li ki + ' ' ' + /uj_|_ r _l, A^_|_ r + • ■ ■ + k 



,i2m e 772111 e 



^n:r,r';i = F h + ■ • • fe i+r _ l , fc i+r + ■ • • + fci+r+r'-l , k i+r+r > +■■■ + h- 2 ) , (IV.14b) 

i- 3 ) , (IV.14c) 
, 2;, (IV.14d) 

(IV.14e) 
(IV.14f) 



-f 1 (ki—i, ki + ■ ■ ■ + /ci_|_ r _i, fci_|_ r , fcj_)_ r -)-i + ••• + &. 

^n™ = ^ lm (^-3, fci-2, h-lih + + • • • + fe,_ 4 ) 



corresponding to the kinematics depicted in fig. 15. The labels on the momenta are defined mod 
n. 



\ i+r+r 1 



i+r 



p4m 

(a) 



i+r+r +r 



i+r i+r 



i-l 



(b) 



i-2 



' (c) 



(-1 



.i+r 



k. 



,2me 



. 1-1 



(d) 



i-3 



.1-2 



(e) 



1-1 



T^Om 



(f) 



Figure 15. The kinematics of the box functions defined in eq. (IV. 14). 



For the purposes of this paper it is useful to exhibit the poles in eq. (IV.8) contained in the 
integral functions with n-point kinematics, 



I 



4m 

4:r,r / ,r // ,i 



r3m 

4:r,r',! 



j2uih 



'4:r;i 



r2me 
L 4:r:i 
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-2r T 



-2r T 
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9 n:i 

L i-3 L i-2 



(IV.15a) 
(IV.15b) 
(IV.15c) 
(IV.15d) 
(IV.15e) 
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where the dimensionful prefactors have been extracted. 
IV. 4 Higher dimension scalar integrals 

Due to the modified integral reduction procedure discussed in appendix I, higher-dimensional 
bubble, triangle and box functions may also appear in the basis of functions in terms of which 
amplitudes are expressed. These integrals may be obtained from the D = 4 — 2e integrals using the 
recursion formulas 

1 / n 

T D=6-2e _ 1 0T _V^. r (i) 



(2/ n - f> x ) , (IV.16) 



(n - 5 + 2e) c 

rD=S-2e _ 1 f 9^=6-26 _ ,-(t),r>=6-2 e \ fIy 1? x 

7 " " (n-7 + 2e) Co V " h ' [ ' 



i=l 
n 



where 



and 



n 



c = J> = E Sy 1 , (IV.18) 

j' = l i=l i,J = l 



Si, = - pj-J 2 , (IV.19) 

is a symmetric matrix. As defined previously the (n— l)-point integral /^Ij is obtained by removing 
the internal propagator between external lines (i— 1) and z from the original n-point (scalar) integral. 
The D = 6 — 2e bubble function is of particular interest and is 

I?--^)=- H1 _Z ){3 _ 2e) (-Ky-'. (IV.20) 

For one- and two-external mass triangles eq. (IV.16) is ill-defined because the Sij matrix is not 
invertible; in this case the D = 6 — 2e integrals may be obtained by direct integration [23] yielding 

,,'-*«■(*?) _ 2e(1 _^ (1 _ 2e) (-Kfr , (W.21a) 

Observe that all higher dimension scalar integrals are related to D = 4 — 2e dimensional scalar 
integrals by coefficients which depend on e. For the integrals which may be obtained from the 
recursion formulas (IV.16) and (IV. 17), the e-dependence appears in the overall factor. Note also 
that the six-dimension triangles (IV. 21) may be expressed in terms of the four-dimensional bubble 
function (IV. 2) with rational coefficients containing explicit e-dependence. 

As discussed in appendix I all e-dependence in integral reduction coefficients may be eliminated 
by including the D = 6 — 2e triangle functions and D = 8 — 2e box functions in the basis of integral 
functions in which amplitudes are expressed; the e-dependence is absorbed into the enlarged set of 
integral functions. Observe that the D = 6 — 2e scalar box is both infrared and ultraviolet finite we 
may drop the e in the overall coefficient in eq. (IV.16) with no effect through O(e ); thus /f )=6_2e 
may be written in terms of D = 4 — 2e integrals with no e appearing in the coefficients. There is 
thus no need to include if =6_2e in the basis of integral functions. 
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